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: CALIBRATION OF A SUBMERGED BROAD CRESTED WEIR 


Richard A. Smith,! M. ASCE 


SYNOPSIS 


work is different from previous publications on the problem in that 

3 the calibration, or correlation of variables, even though the weir is 
ely submerged and tailwater phenomena control the discharge. Pre- 
ritings treat the calibration for the condition of free over-fall, but for 
merged condition, the discharge coefficients were shown to vary over 
range in unpredictable fashion. This paper shows that for free over- 
discharge is a function of the head on the weir, but for the submerged 
m the discharge is a function of the degree of submergence. After 
ion, the computed discharge, when compared with the observed dis- 
showed a probable deviation of approximately 2 percent. 


INTRODUCTION 


es of floods in Louisiana entail the computation of discharge over 
d levees or roadbeds. Available formulas appeared inadequate for 
ng the discharge over the levee or roadbed when considered as a 
rested weir. With available weir formulas, their accuracy was vitally 
nt on empirical coefficients that were little understood in many in- 

In application, a value was somewhat arbitrarily assumed without 
whether the chosen value was even approximately correct. 


at of the Problem 


2 request of the State of Louisiana, Department of Public Works, hy- 
model studies were instituted. The purpose was to calibrate the 
d roadbed, acting as a weir, for two conditions of flow: first, when 


scussion open until August 1, 1959. To extend the closing date one month, a 
n request must be filed with the Executive Secretary, ASCE. Paper 1960 is 
f the copyrighted Journal of the Hydraulics Division, Proceedings of the 
ican Society of Civil Engineers, Vol. 85, No. HY 3, March, 1959. 

and Acting Head, Dept. of Civ. Eng., Louisiana Polytechnic Inst., 


n, La. : 
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the weir controlled the discharge; and second, when the tail-water was 
so that control of the discharge had passed to conditions downstream fr 
the weir. It was assumed that in prototype studies, records of floods w 
yield elevations of the weir crest, head-water, and tail-water surface a 
vealed by high-water marks. 


Objectives in Analysis 
In the analysis of observed data, the specific objectives were threef 


a) Determine a criterion by which it could be ascertained when cont 
the regimen passed from the submerged weir to tail-water condit 
downstream. 


b) Calibrate the head-discharge relationship for the condition when 
trol of the regimen was at the weir. 


c) Correlate the variables, head-water elevation, tail-water elevati 
depth of flow over the weir crest, and rate of discharge, for the ¢ 
tion when control of the regimen was downstream from the weir. 


Description of Models 


Previous observation had shown that in submerged broad-crested wi 
trapezoidal in cross section, the weir coefficient varies as a function o 
slope of the upstream and downstream faces. Therefore, model studie: 
conducted for two models having different slopes. One had slopes of 1 
the other had slopes of 1 vertical on 2 horizontal, with the upstream sl 
being the same as the downstream slope in either case. Other features 
mon to both models, were as follows: Length of weir normal to the dir 
of flow, 4.95 feet; width of weir crest parallel to the direction of flow, 
feet; height of weir crest above the bottom of the channel, 0.425 feet. 
roughness was established by sprinkling the models with masons sand | 
after they had been painted. A diagram of the test apparatus is shown. 
Figure 1. 


Model Operating Procedure 


Head-water and tail-water elevations were controlled by a supply vi 
and tail gate, respectively. Head-water and tail-water elevations were 
served by means of piezometric connections, 1/2 inch diameter, to stil 
wells equipped with hook gauges. The depth of flow over the weir cres 
observed by means of a point gauge located slightly downstream from t 
centerline. This point gauge was read when the fluctuating water surfz 
was above the point and below the point for approximately equal interv: 
time. (Piezometric connections were not deemed feasible at the weir ¢ 
because of the relatively high velocities encountered which could entail 
roneous readings due to stagnation pressures, eddies, etc.) Discharge 
observed by means of a weighing tank and electric stop-clock, which r 
directly to one-tenth of a second. The model was operated throughout 
range of discharge capacity available, up to 1.0 cfs on the weir five fee 
or 0.2 cfs per foot of weir length. At each rate of discharge, the tailw: 
elevation was varied through a full cycle of progressively higher and p 
sively lower elevations in successive runs. 
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s and Nomenclature 
Total discharge rate 
Discharge per foot length of weir 
Coefficient of discharge 
Length of weir, normal to direction of flow 
‘Total head, above the weir crest 
Depth of flow over the weir crest 
Tail-water elevation above the weir crest 
Gravity constant 


Mean velocity 
v 


Froude number, /gy— 
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= Degree of submergence of weir 
ME 


H Relative depth at the weir crest 
Observed Data 


Tables 1 and 2 list data observed for the models with slopes of the 
faces 1 on 1 and 1 on 2 respectively. 


Test Results 


Regimen Criteria 


The criteria were noted from the observations listed below for det 
ing the circumstances when control of the regimen passed from the si 
merged weir to tail-water conditions, or conversely: 


a) When the weir was controlling the regimen, the tail-water eleva 
could be changed without altering the depth of flow over the weil 
or the head-water elevation, with the discharge remaining const 


b) When control of the regimen was in the tail-water conditions do 
stream from the weir, each change in tail-water elevation prod 
corresponding change in the depth of flow over the weir crest al 
head-water elevation, with the discharge remaining constant. 


c 


~~ 


Between successive runs, when control passed from the weir to 
water conditions, or conversely, this change in regimen was re: 
by marked changes in the Froude number of flow over the weir 


d) Data from either model indicated that the shift in control from ' 
to the tailwater, and conversely, occurred in the region of Frou 
bers approximating 0.78 for flow over the weir crest. The high 
bers correspond to the control at the weir, while numbers less 
0.78 correspond to control in the tail-water conditions. 


Discharge Coefficients, Weir Controlling 


From Figure 2, the head-discharge relationships for the models w 
bined with the weir formula: 


3) 
(0) (Cis BS 


Writing the resulting equations with the discharge coefficient as the 
dent variable, this yielded the following relationships: 


Ces Sa ee 
for the model having 1 on 1 slopes and, 
Cag 4” 


< 
for the model having slopes of 1 on 2. A graph of these coefficients i 
in Figure 4. 


TABLE 1 


OBSERVED DATA FROM DISCHARGB OVER MODEL OF ROADBED 


TOTAL 
HEAD 
H 


0.034 
0.033 
0.053 
0.053 
0.053 
0.053 
0.059 
0.123 
0.192 
0.068 
0.052 
0.074 
0.074 
0.072 
0.075 
0.104 
0,073 
0.072 
0,089 
0.089 
0.089 


0.091 


CRESTED WEIRS 


HAVING SIDE SLOPES OF 1 ON 1 


CREST 
DEPTH 


wey 
0.025 


0.025 
0.039 
0.038 
0.039 
0,040 
0.053 
0,123 
0,192 
0.064 
0.038 
0.049 
0.050 
0,050 
0.063 
0.101 
0,052 
0.050 
0.065 
0,065 
0,065 


0.068 


TAILWATER DISCHARGE 


ELEVATION PSR FOOT COEFFICIENT 


Zz 


0,000 
0.000 
0.000 
0.000 
0.023 
0.033 
0.050 
0,121 
0.191 
0,061 
0,000 
0.000 
0.000 
0.009 
0.057 
0.096 
0.039 
0,000 
0.000 
0.000 
0,000 


0,061 


ee dan eet 
0183 


0180 
0369 
0364 
0358 
0356 
0354 
0351 
20351 
.0352 
0351 
0610 
0609 
0610 
0608 
0611 
0610 
0609 
0840 
0841 
0841 


. 0838 


OBSERVED 


C 


2.92 
3.00 
3.02 
2.98 
eK) 
2.92 
2.46 * 
0.815 * 
0.417 * 
(1.99)* 
2.95 


3.04 


3.04 


TABLE 1 (CONTINUED) 


RUN 
NO, 


23 
24 
25 
26 
27 

28 
29 
30 
on 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 


45 


CREST 
DEPTH 


Is 
0,084 


0,067 
0.066 
0.074 
0.074 
0.074 
0.074 
0,096 
0.125 
0.125 
0,082 
0,080 
0,079 
0,078 
0.076 
0.074 
0.073 
0.084 
0.084 
0.118 
0.154 
0,236 


0.112 
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TAILWATER DISCHARGE OBSERVED 


BLEVATION PER FOOT COBFFICIENT 


z 


0,077 
0.052 
0,044 
0.000 
0,000 
0.000 
0.000 
0.087 
0.118 
0.118 
0.078 
0,070 
0.062 
0.051 
0,025 
0.010 
0,000 
0.000 
0,000 
0,008 
0.148 
0.233 


0.100 


eee) Fae 
0839 


0842 
0839 
- 1006 
- 1004 
- 1006 
- 1005 
-1002 
0976 
- 1003 
- 1006 
-1011 
-1011 
- 1009 
- 1009 
1011 
- 1008 
e1411 
- 1400 
-1411 
1411 
- 1404 


-1407 


C 


2.82* 
3.12 
3.16 
3.18 
3.18 
3.18 
3.18 
2.79% 
.2.08* 
2.17* 
3.00* 
3.05* 
< wt 
3.19 
3.19 
3.20 
3.19 
3.16 
3413 
2.78* 
2,19* 
(1.22)* 


2.84* 


HY 
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TABLE 1 (CONTINUED) 


) TOTAL CREST TAILWATER DISCHARGE OBSERVED PROUDE 
| RUN HEAD DEPTH BLEVATION PER FOOT COBFFICIZNT NUMBER 
| NO. H 
| awe Bo Mas: SERRE SSS a es Wee Cc eb eg 
| 46 0.129 0.107 0.089 £1417 3,06* 0.714 
47 0.126 0.088 0.054 .1405 3.14 0.950 
48 0.126 0.084 0.019 .1409 3.15 1.021 
49 0.126 0.084 0.003 .1410 3.15 1.022 
50 0.126 0.084 0.000 .1410 g205 1.022 
51 0.126 0.084 0.000 1426 3.19 1.033 
52 0.148 0,098 0.000 .1832 3f2i 1.053 
53 0.149 0.097 0.000 1831 3.19 1.070 


. ) Observed coefficient is markedly different from computed mean, but 
the observation is in the range of small heads where large percent 
errors may be anticipated. 


* The weir is drowned out, and the head-discharge relationship is 
controlled by tail-water phenomena. 


charge Coefficients, Tail-Water Controlling 


n Figure 3 is shown graphs of the degree of submergence and relative 
th at the weir crest as related to the Froude number. The equations of 
curves and the graphs were derived by the method of least squares. In 
er to avoid undue confusion of the graph, the individual data points are 
tted. 


Model with Slopes of 1 on 1 


‘or the model having slopes on the weir faces of 1 on 1, the degree of 
mergence is related to the Froude number by the equation, obtained by 


tt squares analysis: 


(2) + 4 = 980 (1) 
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TABLE 2 OBSERVED DATA FROM DISCHARGE OVER MODEL OF ROADBBD 
HAVING SIDE SLOPES OF 1 ON 2 


TOTAL CREST TAILWATER DISCHARGE OBSERVED 


RUN HEAD DEPTH ELEVATION PER FOOT COEFFICIENT 
NO. H PSs See eo eee eee 
61 0.039 0.024 0,000 0218 2.83 
62 0.035 0.021 0,000 .0188 2.89 
63 0.053 0.033 0.000 .0349 2.86 
64 0.052 0.033 0.000 0345 2.90 
65 0.052 0.033 0.032 .0337 2.83 
66 0.058 0.048 0.048 0332 (2.37)* 
67 0.069 0.045 0.000 .0527 2.91 
68 0.069 0.049 0.046 0528 2.92 
69 0.072 0.053 0.055 0529 (2.74)* 
70 0.090 0.058 0.000 .0838 3.10 
71 0.091 0.058 0.028 .0837 3.04 
72 0.101 0.087 0.086 .0837 2.61* 
73 0,101 0.065 0.000 1014 3.16 
74 0.099 0.065 0.000 .0987 3.16 
75 0.099 0.063 0.005 0989 5.37, 
76 0.099 0.064 0.041 .0987 3.16 
77 0.119 0.108 0.106 0989 2.41% 
78 0.099 0.065 0.055 .0990 3.17 
79 0.126 0.082 0.015 .1407 3.14 
80 0.126 0.082 0.066 1412 3.16 
81 0.149 0.136 0.134 1406 (2.44)* 
82 0.141 0.127 0.123 1410 2.66 * 


TABLE 2 (CONTINUED) 


RUN 
NO. 


83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
100 
101 
102 
103 


104 


(105 


CRESTED WEIRS 


TAILWATSR DISCHARGE 


ELEVATION PER FOOT COEFFICIENT 


Zz 


0.104 
0.089 
0.069 
0,000 
0,093 
0,070 
0.022 
0.121 
0,091 
0.079 
0.059 
0.046 
0,000 
0.081 
0,064 
0.043 
0.000 
0.091 
0.060 
0,028 
0,090 
0.065 


0.011 


oes le 
-1410 


-1409 
1411 
- 1410 
21822 
1781 
-1822 
. 1008 
. 1007 
- 1008 
- 1009 
- 1008 
- 1009 
0858 
-0843 
0803 
20804 
0631 
~0622 
.0620 
0421 
-0420 


20395 


OBSERVED 


Cc 


2.97* 
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TABLE 2 (CONTINUED) 


TOTAL CREST TAILWATER DISCHARGE OBSERVED FROUDE 
RUN HEAD DEPTH ELEVATION PER FOOT COEFFICIENT NUMBER 
NO, H y. z g C F 
106 0,053 0.036 0.001 .0395 3.24 1.019 


( ) Observed coefficient is markedly different from computed mean, but 
the observation is in the range of small heads where large percent 
errors may be anticipated, 


* The weir is drowned out, and the head-discharge relationship is 
controlled by tail-water phenomena. 


The relative depth at the weir crest is related to the Froude number by 
the equation, obtained by least squares analysis: 


2 
y eee 
ee PES OA ee eS OOS ( 
H ay. 


Combining equations (1) and (2) with the continuity equation, 


C Fey. 


the resulting equation is as follows: 
3 2.2 : 2.2] 3 | 
g-if[g[90-(2)"]lecsvsz(Z)P 
By comparing this equation with the weir formula, 


Ochs 


it is seen that, 
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Ce | 9 9e-(Z)"] [6s risen 5,27) ars (4) 


Model with Slopes 1 on 2 


‘or the model having slopes on the weir faces of 1 on 2, the degree of sub- 
gence is related to the Froude number by the equation, obtained by least 


res analysis: 


ze ae 
aH + Baga = .965 (5) 
A Cae 


N|z 
tae 3 
Soh 
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Sa | 
BM 
We ae Sh Na 
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The genes ated Bie des Sete apni ude number t 
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2 
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os 


When equations (5) and (6) are combined with the continuity equation, 
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“pag ane 


; resulting equation is as follows: 


. 2 7 oe eee. 3 
| g=H2/9 2.a6- 2.962] [027 x 992 | (7) 


comparing this equation with the weir formula: 


@—-CLAe 


, 


Ss seen that, 


: 3 
Cx [yg] 2.86 - 2.96 | [027 +9912] (8) 


Thus it is seen in equations (4) and (8) that the coefficients of discharge 
either model is a function of the degree of submergence (z/H). It should 
noted that this correlation is for the condition of the weir drowned out 

| control of the regimen is in the tail-water. A graph of these coefficients 
shown in Figure 4, representing equations (4) and (8), above. 


servations Discarded as Erroneous 


In comparing observed data with the computed mean values, it was noted 

t large deviations occurred in 2 runs out of a total of 53 for one model and 
! runs out of a total of 46 for the other model. In each instance of large 
iation, the total head from head-water to tail-water was 0.018 feet or less. 
ure 5 shows that for smaller heads, an error in observed head will pro- 

e a larger percent error in computed discharge. Therefore, in arriving 
he probable deviation, the larger deviations were omitted from the 

lysis. Those runs omitted are indicated in Tables 1 and 2, by parentheses 
losing the numerical value of the coefficients in Figure 4. 


ors in Discharge Relative to Errors in Head 


The following argument will show the relative deviation in computed dis- 
rge caused by a corresponding deviation in observed head. The weir 
mula, 


3 
Qieipcye Ashe 
7 be expressed logarithmically, 


ha) Ce FRG. ee In L + #InH 


and L are treated as constants then the derivative of this equation will 
d, 


3 
CLF 


COEFFICIENT “C’' IN @Q 


ate 9 8 7 6 O5 


March, 1959 Hy 


HEAD-ON"“WEIRIN FEE a 
WEIR CON TROLLING 
0.0 s/ iz a 4 ew 


£ 


( 


DEGREE OF SUBMERGENCE (% ) 
TAIL-WATER CONTROLLING 


FIG. 4 WEIR COEFFICIENTS 
FOR MODELS OF ROADBEDS 
WITH CONTROL EITHER AT THE 
WEIR, OR IN TAIL-WATER 
PHENOMENA. 
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| small, finite increments, this may be expressed as follows: 
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‘The percent deviation in computed discharge is 1.5 times the percent 
yiation in observed head. This relationship is shown graphically in 
sure 5. 
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Deviation of Observed Data From Computed Mean 4 

The laws or probability were applied to compare the observed discharge 
with mean values computed by the weir formula and using coefficients from 
Figure 4. This comparison indicated that the observed data agreed with th 
computed mean values with a probable deviation of a plus or minus 1.4 per. 
cent for the model with the slopes of 1 on 1. Similarly the probable devia- 
tion was a plus or minus 2.1 percent for the model with slopes of 1 on 2. 
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TABULAR SOLUTION OF OPEN CHANNEL FLOW EQUATIONS 


Henry A. Babcock,! A.M. ASCE 


ABSTRACT 


This paper presents and explains a table by which cubic equations which 
ise in open channel flow problems may be solved without the tiresome trial 
d error method usually necessary. 

Use of the table is limited to the special case of rectangular channels. 


Most engineers will cheerfully admit their inability to solve equations be- 
nd quadratics by any method other than “cut-and-try.” Special methods for 
lving cubic equations have been devised for particular equations of frequent 
currence, and such a method is presented here for an equation that arises 
open channel hydraulics. 

The familiar specific-energy diagram forms the basis for the solution. 
though this diagram is reproduced in practically every fluid mechanics 
tbook on the market,(1) none of them explain the use of the diagram in 
actical calculations. Figure 1 shows a specific energy diagram plotted to 
ule in dimensionless form. Actually, a diagram cannot be read with suffi- 
nt precision for many purposes and so a table giving the same data is also 
sented. 4 

A brief summary of the theory follows. 

For a rectangular channel, the specific energy at any cross section is 


Eis PP RAS (1) 
2 


i the critical depth is 


te: Discussion open until August 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1961 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
“American Society of Civil Engineers, Vol. 85, No. HY 3, March, 1959. 


Asst. Prof. of Civ. Eng., Colorado School of Mines, Golden, Colo. 
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where the symbols have the following meanings: 


Eg = specific energy 

d = depth of flow 

q_ = rate of flow, in cfs per foot width 
g = acceleration of gravity 

d. = critical depth 


Division of Eq. 1 by Eq. 2 gives 
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E 
Hpiues of = for values of a € 4.99 are given in Table 1. 
c Cc 


Figure 1 was plotted from the data in Table 1. 
Use of the diagram will usually yield results of satisfactory accuracy, 
wever, greater mathematical precision is obtainable from the table. 
Problems are solved by writing Bernoulli’s Equation in the usual way, and 
n by appropriate manipulation reducing it to the form of Eq. 3. When this 
one, one side or the other of Eq. 3 can be evaluated from given data, and 
n the other side can be found in Table 1. Finally, multiplication of tabular 
lues by de will yield a numerical result. 
If it is necessary to include an energy loss between two cross sections, 
table or graph may be used by making a slight modification of the theory 
lined above. Say that between the upstream section and the downstream 


Vv 
ction there is a loss equalto K ae . (The subscript 2 refers to the down- 
‘eam section.) Bernoulli’s theorem written between the upstream section 1 


i section 2 is 


dt +k (4) 


(a 
q 

Suh ae (1s = (5) 
or2 


Setting the derivitive of (5) with respect to dg equal to zero and solving 


dg yields 
3) (i+K) a* 
eae iw er eeeee (6) 
o 9 


This value of dg will be designated d,,. The relationship between dy and 


is 
3] i 
a= (/+ K) a. (7) 


Division of Eq. 5 by Eq. 7 gives 


Fen, 2, 1/2 ve (8) 
A, qd, 2 a> 


This equation is of the same form as Eq. 3 and when either side is known 
other side may be found from Table 1 or Figure 1. 
Illustrative Problems 


The following problems would lead to a cubic equation if the accompanying 
le were not used. Computations are carried to three decimal places to 


ure accuracy of the second decimal place. 
j , 
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TA BL e ft 
E,/d, AS A FUNCTION OF d/d, 


d/de, 0500 © 0.0L ©. 05025250.03%. 205049 0505.) 006" 0-07 0.08 


0.0 5000... 1250.4 5555.6  312.55-200s1 5,138.9) 102,178.21 
Oxl> 50.10. 41.43 634.86 -29.714°25.65" 22537 19,0917 4s eat OL 
O22" 12.70 BU255: 10,5 5009. 082 (0.920 o.200 6 f OIG e iste | 1G. Oe 
0.3 5.856 5.513 5.203 4.921 4.665° 4.432 4.218 4.022 3.843 
0.4 32525. 3.384 35254 3.134 932025) 23919927823) 2.736 ee 3650 
0.56262 25009E2543202. 9699925310025 2552 7203022154 Ze 108 R066 
0.6 1.989 1.954 1.921 1.890 1.861 1.833 1.808 1.784 1.761 
OTcackas2Z0 mL 202. 1 6845 1. 668 01.653 elo 1626 2An Gls bie 
0,8: 1.58, 1.5/2. 1.564: 1. 55Grrb 540 542 Loos Gari ook) eae 
O59" 22507) 1.514 SLE AL 508) US06 1 1 504. 61.502 1 0507 
L260, 11350001. 500,” 1.50) <b. 50 APSO, 5h2504° 91 5050: t. 50). oue 
Le PL SES) ASS YG: St S5L on S22 aco E526 S153 2e aS oes Ooo 
U2 E5547. 552 “15556 21, 560501 .565™.1.5/0_ 1.575. 1 S80 Rain seo 
1230. 1.596- E.60L Ua 607523613, 2 618 15624-1630 1636 aeb47 
4 12655 ‘L662 1.668 51674 16st 1.688) 1.695 st fOr a 270s 
TS Le 22 wa 2921796 2 7a LOSE LP oo, ok OO ek to oO 
156°). 1.795..4.803. Ps810.22.818 1.826’ 1.834. L841 1.849 91857 
Ld ebaOts. Le 80l el, 889 | Loot a) 90s e9ls)) J. 92 1 Osorno sae 
Ps8> 1.954. 12963 LOPE O7 oe. Le S88rs 996) 22 04ers Ol siaes02Z 
iD) 5 2.038% 2,047 (27056: "2.06402 .073: 2.062 (2.090). 2.099" 25808 
20. | 2Zel2Z5. -25134~ 2.1422 2.151 ene. 160) 26169 27S 2 Lop aieed 96 
2m 2.21375 25222 2.23.2 261 2249. 22258 82 26722708 ees 
fee Sep AOS) 2adi se 2 522 B2.590) 2a SW 2eS49 es ee 550 Le SOT mee 
2.3» 2394. 2,404 2.413" 2.422 2,431 2,.440..2.450 9.2.4595 2-465 
Dak GestOh  SetI0”  ee509 per Oka, An dek 26550 Bee Das eee al OL 
B50 2.580) 92.589" 2.599" 25608! 2618 25627 723636) 2.646) 922655 
20) 2074 22083 25099 2510. 2 Foe TelL nN Cel ol eee PU0me2. FOO 
Zed) ie OD es Lad Sar ei OOD Lote ae OO EaeisOl Ole CLO fel OI Ome Boas 
Ce Oe 26004, 2.073 2,003) e.00e) 2902 ee. Ile 2 92) 62, 93b 2 ORO 
BeBe 26.900) 1229699 52.9795 22988) 27998 3008  S20L7 FS027 5 220s6 
37091 360560.3 065.035 070523.084 3.09453 L043 VIS 223 2s so aS 
Sele peo oe so LOLs Souk Selo 25,19 e200) os. 2) Oeeae2 20naia 229 
Bie e249 poe eos oe eOR | Se2sB eo 2OGN sco fea. S07 Oh Tus 2G 
Sede 5 Sad4Oe 3e300™ Oo. 905 23.37 DwowS0De 3 uO9D™ e404 ea ee wes ed 
324 3.443) 3.453: 3.463" 93. 492093 482 23.402 — 350283 512 “37520 
SoD) 735.941) 355506 53. 560 355708" 3\560e 3. 590 > 3600) 63). 6094 13.619 
36) “36395 (8.1648) 35658) 34.668, 635678) 3. 6882-3.697 232707 easeva 7 
Sil SST Ve VEOr Se /D0 Se tOO 3.700 A100 aoe OS ms MOUS mn SOL 
3.8 3.835 3.844 3.854 3.864 3.874 3.884 3.894 3.903 3.913 
329 35933; °3.943% 3,953. <3. 962) 35972) -3.982 . 3992) 4.002) 4 Oke 


4.031 4.041 4.051 4.061 4.071 4.080 4.090 4.100 4.110 
4.130 4.140 4.150 4.159 4.169 4.179 4.189 4.199 4.209 
4.228 4.238 4.248 4.258 4.268 4.278 4.288 4.297 4.307 
4.327 4.337 4.347 4.357 4.366 4.376 4.386 4.396 4.406 
4.426 4.436 4.446 4.456 4.465 4.475 4.485 4.495 4.505 
4.525 4.535 4.544 4.554 4.564 4.574 4.584 4.594 4.604 
4.624 4.634 4.643 4.653 4.663 4.673 4.683 4.693 4.703 
4.723 4.732 4.742 4.752 4.762 4.772 4.782 4.792 4.802 
4.822 4,832 4,842 4.851 4.861 4.871 4.881 4.891 4.901 
4,921 4,931 4.941 4.951 4.960 4.970 4.980 4.990 5.000 
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ample 1 


_At what two depths could a flow of 400 cubic feet per second at a specific 
ad of 7 feet be carried by a rectangular channel 9 feet wide ?(2) 
From the given data compute the following: 


: 3/44 4aé 
: er pia as 9 894 Suh t 


= 1.774. 


E 
From Table 1, for this value of —* the two values of a. are a = 1.571 
de de de 


d 
d do = 0.674. Multiply these values by d, to obtain 
d, = 6.198 feet (tranquil flow) 
d, = 2.659 feet (rapid flow) 


cample 2 


A channel 10 feet wide and 6 feet deep carries 300 cubic feet per second. 
transition reduces the width to 9 feet with no change in bottom elevation. 
nd the depth in the narrow section. 

From the given data compute the following: 

In the 10 ft. wide section 


! rz) 
3O 
a, = \/S =3.036 ft. 
g 
Ope pee Se = (976 
a, 3.036 
; dy oe Eaol 
Enter Table 1 with this value of a. and obtain rope = 2.104. Hence 
c c- 


-1 = 2.104 x 3.036 = 6.388 feet. 
In the narrower section 


2 
pm LOSS |g 305 7- 
C-2 g 


Es-2 6.388 
If losses are neglected Eg_9 = Eg-j tether? and eae = 3-957 = 1.961. 


om Table 1 the corresponding value of Tan 1.808 from which 
| c 
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dg = 1.808 x 3.257 = 5.889 feet. 


Example 3 


If a rectangular channel 8 feet wide carries a flow of 250 cubic feet per 
second at a depth of 5 feet, what change in surface élevation will be produce 
by a local rise in floor level of 1/2 foot? 


3 2 
d= ee ©3118 FF. 
9 


See Se ee EL 


d E 
Enter Table 1 with this value of = to obtain a ="1:797. "From 
c 
Bernoulli’s theorem S 


Es-1 = Eg-2 + 0.5 
Divide both sides of this equation by d, and rearrange the terms to obta 


All quantities on the right of this equation are known. 


ar = £79 Zaz, 160 * 1263.0 
a, 


d 
Enter Table 1 with this value and obtain + = 1.371 from which dg = 
c 


4.28 ft. Hence the surface will drop 0.22 ft. 
Exanple 4 

The bottom of a rectangular fli me is 10.00 feet below the pool level in a 
reservoir. If the entrance loss is 0.2 A find the water depth in the flume 
for a flow of 70 cubic feet per second .. foot. 

Write Bernoulli’s theorem between the pool surface and the channel 


2 2 
V V 
Me ® d* 55 * 0255 
Substitute V = 4 and combine the last two terms. 


2 
H= de La > 
2gd 
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Sf-ri sat 
md. = as and divide each side by dy, The result 


| Hd » £ (2a) 
| poe eS bd 


This is the same form as Equation 3 and Table 1 may be used with din 
lacing d.. Therefore, compute 


3 2 
70 

ON ihe) mma S ‘ 

Ln, Sica 5.669 ft 

H i /0 one 

por Reo 4 
Enter the table with this value of a and obtain a: = 1.558 or 

aa dm 
2.558 x 5.669 = 8.83 ft. 
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RELATIONSHIPS BETWEEN PIPE RESISTANCE FORMULAS 


Walter L. Moore,! A.M. ASCE 


SYNOPSIS 


The relationship between modern concepts of pipe resistance and the older 
pirical formulas is clarified, and a simple procedure developed to derive 
exponential formula applicable to a known range of flow conditions. The 
alysis makes apparent the range of flow conditions for which any particu- 
‘exponential formula is applicable. Limitations of the equivalent pipe con- 
tt are discussed. Data on resistance measurements in water mains in 
rvice indicate that the head loss varies with the discharge to a power near- 
equal to 2.0 rather than 1.85 as commonly assumed. 


INTRODUCTION 


Methods for predicting the resistance encountered by a fluid flowing 

ough a pipe have long been of concern to mathematicians and engineers. 

= problem has been recognized as an important one in the design of water 
ply lines and water distribution systems as well as in process plants of 
‘ious types. Recently the use of long pipelines for transporting large 
umes of fluids has renewed interest in problems of pipe resistance. 

Over the years, different approaches have been made to the problem of 
dicting pipe resistance and a multitude of equations and graphs have been 
ised for evaluating the resistance. In recent years, the literature on the 
ject presents semi-rational approaches leading to relations applicable 

Tr a wide range of conditions. These approaches have been offered to re- 
ce the older empirical relations. Many practicing engineers, however, 

i to the old familiar relationships without a clear realization of their 
itations. Perhaps this should not be surprising because so little has been 
tten to explain the interrelationship between the various methods of 


8: Discussion open until August 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1962 is 
gart of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 85, No. HY 3, March, 1959. 


Prof. of Civ. Eng., The Univ. of Texas, Austin, Tex. 
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evaluating resistance. This paper attempts to clarify this interrelationshi 
showing exactly how the older empirical equations may be regarded as spe 
cial cases of the more recent general resistance relations. It is hoped that 
with clarification of this interrelationship more use will be made of the re- 
cent information as it is presented in the general resistance diagram. This 
information may be used in equations of the same exponential form as the 
common empirical formulas by means of the principles presented herein. 


Development of Pipe Resistance Relations 


A brief review of some selected major steps in the development of pipe 
resistance relations will help to clarify the resistance relations to which 
reference is made. In Table I, some selected major steps in the develop- 
ment of pipe resistance relations are shown giving the approximate date of 
the development, the name of the investigator commonly associated with it, 
the resistance relation, and the conditions for its application. 

Some of the resistance relations, namely the Chezy, Manning, and Hazen 
Williams equations, are customarily in a form convenient for solving for th 
velocity V. In Table I, these are also written in a form in which the head 
loss hg is expressed in terms of the other variables after substituting for t 
slope S as o . In this form, they are more comparable to the common fo! 
of the Darcy-Weisbach equation. The relations listed as previous to 1930-: 
in Table I are all based on empirical studies with the flow of water. They 
can all be expressed in the general form 


hy = CL ym ( 


R 2 


where C is a coefficient to account for the pipe material, L the length of t 
pipe, R the hydraulic radius, V the mean velocity, and z and m are expo 
nents. More generally they could all be expressed as 


h, = KQ ( 


where K is a coefficient which depends upon the length, diameter, and ma- 
terial of the pipe and Q is the discharge. Equations of this type are terme 
exponential equations because the head loss is expressed in terms of the pi 
characteristics as embodied in K and the discharge Q raised to some expo- 
nent m. 

The modern understanding of pipe resistance began about 1930 with the 
Prandtl-Karman relations which were supported by the experiments of 
Nikuradse, Colebrook, and White. These modern relations are based on an 
understanding of the variation of the resistance coefficient f in the Darcy- 
Weisbach equation. For the steady uniform flow of an incompressible fluic 
the complete picture is well summarized in a general resistance diagram_ 
such as proposed by Rouse and by Moody in slightly different forms. The 


values of f from the general resistance diagram are substituted in the Daz 
Weisbach equation 


TABLE I 
Selected Major Steps in the Development of Pipe Resistance Relations 
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Sn =e = erie 
Year Name Relation Conditions for Application 
= iF 
1775 Chezy V =CNVRS Flow of water in open channels. 
= Lv2/rce2 
h, = LV /RC 
Seem 
1€45 Darcy= ee een v? All pipe flow. 
Weisbach fy oD oe 
1859 Poiseuille h, ‘. ae Laminar flow in pipes. 
1883 Reynolds Criteria for Laminar Flow in pipes. 


vD 
Ny = 37 < 2000 


1869 Manning y - Lhd 72/3512 Turbulent flow of water in 
n channels and pipes. 


1892 Freeman Exhaustive tests and Flow of water in commercial pipes. 
tabular results. 


114 Hazen- 


0. 63,0 “54 
Williams 


Ve 1.318C,R Turbulent flow of water in pipes. 
+85 


Hel 27,809,115 
=e ea ee 


1930-31 Prandtl- Turbulent flow of any fluid in 


a 
Karman = 2 log N, VF - 0.8 a "smooth" pipe." 
Te bP Shite BP 
1 To 
to 2 log + 1.74 Turbulent flow of any fluid in 
f a "rough" pipe. 
1932-33 Nikuradse - Systematic test with Pipes artifically roughened with 
controlled roughness. sand grains. 
1939 Colebrook- at ae 2.51 e Commercial pipe, turbulent flow, 
White mare (Sh aint oF "smooth," "rough," and transition 
Je ™ cS chalph conditions. 
19he Rouse Plot of f = General resistance diagram for 
(at Jz,£ ) commercial pipe: Laminar and tur- 
(ly Np Vf bulent flow, smooth, rough, and 
transition conditions. 
ES SSS Se 
19h} Moody Plot of f = "Moody Diagram” - for commercial 
e pipe: Laminar and turbulent flow, 
# 95 ) smooth, rough, and transition 
conditions. 


f Deze 


raluate the flow resistance. This relation may also be expressed as 


Z (4) 


e K, depends on the length, diameter, and roughness of the pipe as well 


ie Reynolds number of the flow. 
ach of the two types of formulas, the exponential type and the one based 
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on the general resistance diagram, has some advantages as well as disad- 
vantages. With the exponential type, the calculations are convenient; and | 
this reason it is widely used by many practicing engineers. The calculatic 
are direct, avoiding the need for trial-and-error solutions. The reliabilit 
of the calculated resistance is questionable, however, because of the unce1 
tainty in selecting the applicable equation from the very large number of e 
ponential formulas that have been proposed. As late as 1951, Blair(10) rej 
to the multitude of formulas already proposed and proceeds to suggest fou: 
new exponential formulas designed to cover a wide range of pipe materials 
Again in 1954(14) these formulas were reexamined and eleven new exponer 
formulas (divided into four groups) were proposed to represent a wide ran 
of conditions. Because the exponential formulas contain no term involving 
the fluid viscosity, they are necessarily limited to conditions where viscos 
has little effect or to a restricted range of Reynolds numbers. 

The modern general resistance diagram, however, accounts in a ration; 
manner for all of the factors which influence pipe resistance. A clear und 
standing of the mechanism of pipe resistance has come from this approacl 
and the general resistance diagram clearly shows the effect of the various 
factors entering into the problem. The improved knowledge incorporated 
the general resistance diagram, however, has not received as wide usage 
among practicing engineers as it deserves. This is probably largely duet 
its inconvenience as compared with the simpler exponential formulas. Wi 
the general resistance diagram in its common form, trial-and-error solut 
are required for finding the rate of flow or the size of pipe. 

Because of the convenience of the exponential formulas and their wide- 
spread usage, it is desirable that the interrelationship between them and t 
general resistance diagram be clearly understood. In this way it may be 
possible to benefit from the convenience of the exponential formulas as we 
as the completeness of the general resistance diagram. The interrelation 
between exponential formulas and the general resistance diagram has not 
previously been presented in complete form. 


Interrelationship Between Moody Diagram 
and Exponential Formulas 


For a limited range of Reynolds numbers and a fixed value of , thep 


of log f versus log Np may be approximated by a straight line and f ma 
be expressed in terms of the Reynolds number Np as 
Cc 
= o 
 (2-m) 
NR 


rh 


where Co is a constant and (2-m) represents the negative slope of the ple 
log f versus log NR on the Moody diagram. If this is substituted into a f 
Darcy-Weisbach equation and the Reynolds number is expressed in terms 


the mean velocity, pipe diameter, and kinematic viscosity of the fluid, the 
results ; 


oon hag Y 


= —-__ sg 
"4 2e(/4)™p 3+m 


m 


Ss 
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— loss may then be expressed as 
m 
| he = KQ (7) 
pre 
| Cy prem) 1 
Ko = —~-——____,____ (8) 
fe) Seat m 


this expression both the coefficient Co and the exponent m will be func- 
ns of the Reynolds number and the relative roughness of the pipe = . Thus, 
en the head loss is expressed in the form of equation (7), the coefficient 
depends upon the pipe characteristics of length, diameter, and relative 
ighness as well as the Reynolds number and fluid viscosity (see equation (8)). 
The exponent m may be conveniently determined from the Moody diagram 
sause (2-m) is the negative slope of a curve on the Moody diagram. Thus 
es of constant m may be drawn on the Moody diagram by connecting points 
constant slope. This has been done in the modified Moody diagram of Fig. 1. 
The modified Moody diagram makes clear why so many exponential formu- 
have been proposed. It is evident that to completely represent the resis- 
ce problem an infinite family of exponential formulas is needed with the 
yonent m varying continuously from a minimum value of about 1.7 for 
ooth pipes at a Reynolds number of about 4000 to a maximum value of 2 

a relatively rough pipe at large Reynolds numbers. For a given value of 
ynolds number and relative roughness the proper value of the exponent may 
read from the modified Moody diagram. 
The use of the proper exponent m is important when the head loss is to be 
ermined for different discharges or when a change of head loss is to be 
ated to a change in discharge as in the Hardy Cross network analysis. Dif- 
entiation of the equation hg = KQ™ to yield ee = me shows that the 
ative change in head loss is m times the relative change in discharge. 
7 of the common exponential equations such as the Hazen-Williams, 
nning, or Darcy-Weisbach (with constant f) may be made to yield the 
rect head loss at a given discharge by proper choice of C, n, or f. How- 
r, the head loss at other discharges as well as the relation between a 
nge in head loss and a change in discharge will be incorrect unless the 
onent m of the selected formula corresponds to the combination of rela- 
roughness and Reynolds number for the pipe under consideration. The 
lified Moody diagram gives an easy check on the proper m value for any 
ibination of relative roughness and Reynolds number. Thus the Hazen- 
liams formula with its exponent of 1.85 corresponds to a combination of 
ative roughness and Reynolds number lying along the line m = 1.85 on the 
dy diagram. The Manning equation or the Darcy-Weisbach (with constant 
orresponds to conditions lying on or to the right of the line m = 2.0. 
With the aid of the modified Moody diagram it is convenient to obtain an 
onential formula for a reasonable range of conditions which will accurately 
ect the resistance characteristics as embodied in the Moody diagram. The 
owing steps have been found most convenient: 
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hh Transition zone 
R= — Complete turbulence, rough pipes = 


7H \ 
PIPE MATERIAL (new) 
brass, lead, glass, transite 
comm. steel, wrought iron 
galvanized iron 
cast iron (avg.) 
concrete 
riveted steel 


Friction Factor 


VD V4R ; £:10 
s—_— . D 
Ne Tl as (ft. sec. units) 


Fig. | Modified Moody Diagram 


. Select approximate values of Reynolds number NR and relative roug 
e 
ness 5. 


. From the modified Moody diagram obtain the exponent m and the fri 
tion factor f. 


. For a velocity consistent with the Reynolds number used in Step 1 ca 
culate the head loss in the pipe using the friction factor f in the equé 


tion. 
hy Spt L We 
D 2g ( 
Compute the discharge Q in the pipe. : j 
h 


f 
4. Solve for Ko F con 
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is gives values for K, and m in the equation hs = Kj)Q™. Because Kp and 
may change with Reynolds number and hence with Q the exponential equa- 
n thus found should be limited to a reasonable range of variation in Q. For 
ditions falling in the transition range between smooth -pipe and rough-pipe 
w on the general resistance diagram, the results from the exponential 
ation agree with those obtained from the resistance diagram within 5 per 
t for a twenty-fold variation in Q. The agreement will be closer for con- 
ions approaching those for rough-pipe flow, that is, m approaching 2, as 

li as for conditions approaching smooth-pipe flow. Examples in the 

pendix show in detail the method of obtaining exponential formulas from 

2 modified Moody diagram. 


Equivalent Pipe Concept 


The equivalent pipe concept is commonly used to analyze the flow for pipes 
Series or in parallel. This may be done by replacing the individual pipes 
the system, for purposes of calculation, with several pipes all of the same 
ameter but whose lengths are chosen to give loss characteristics equivalent 
the original pipes. Two pipes are thus considered equivalent if they have 
ual head loss at equal rates of flow. This may be expressed in terms of 

> exponential equations as follows: 


1. For pipe No. 1 let hg = Ko, @ 1. 
2. For pipe No. 2 let hg = Koo QQ, 


the normal use of the equivalent pipe concept the exponents mj and mg are 
sumed equal when a particular exponential equation is selected for the 
oblem. Then if pipe No. 2 is to be made equivalent to pipe No. 1 the length 
pipe No. 2 is so chosen that Koo is equal to Koy: 


The previous discussion of exponential formulas makes it clear that the 
9onents m may not be equal for all pipes in a system and if this is true 
‘re iS no one pipe which can be made equivalent to the various pipes in the 
stem. Equivalence may be achieved for one rate of discharge, but not for 
road range of discharge. 

The effect of a difference in the exponent m may be illustrated by an ex- 
ple in which a 12-inch diameter steel pipe is compared with a 6-inch 
meter cast iron pipe. A 100-foot length of 6-inch cast iron pipe having a 
sistance relation h¢ = 1.88Q1-95 may be compared with 3660 ft. of 12-inch 
el pipe having a resistance relation hg = 1.99Q1-88, The length of 12-inch 
e was chosen so that the two pipes would be “equivalent,” that is the head 
s for the two pipes would be equal, at a discharge of 2.25 cfs. The results 
this comparison are shown in Table II. It is apparent that while the pipes 
: equivalent at a discharge of 2.25 cfs, at a flow of 0.5 cfs the head loss in 
steel pipe is 12 per cent greater than that in the 6-inch cast iron pipe, 

ile at a flow of 10 cfs the head loss in the 12-inch steel pipe is 11 per cent 
s than that in the 6-inch cast iron pipe. This demonstrates that when the 
onents m differ significantly the equivalent pipe concept may lead to in- 
urate results. : 

Although there are some situations where the variation in the exponent m 
significant, in many applications of resistance formulas for design purposes 
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TABLE II 


The Effect of the Exponent m onthe Equivalence of Pipes 


9 3660 ft. of 100 ft. of 6-inch Ratio 
cfs 


12-inch Steel Pipe Cast Iron Pipe hy Thy 
h, =1.99Q1-88 ft. ny =1.88Q }+ 96, 12 “6 
£ f 
12 6 
(Oras) 0,54 0.49 , ae Bt 
1,0 199 1.88 1.06 
2.0 Key 16:25 1.01 
3,0 15.70 16.00 0.98 
4.0 27.00 28.20 0.96 
6.0 57. 10 62.00 0.93 
8.0 99.50 109.00 0.91 
10.0 151.00 169.00 0.89 


the variation in m will be of little practical significance. If a pipe system 
includes pipes with a broad range of relative roughness and a broad range o 
Reynolds numbers, the conventional “equivalent pipe” concept may not be 
adequate to explain the operation of the system. However, its performance 
could be understood in terms of the applicable resistance equations as dete 
mined from the modified Moody diagram. In many applications of resistanc 
formulas for design purposes, the inaccuracy due to the variation in the ex- 
ponent m will be of little practical significance because the precision of the 
basic design conditions does not justify great refinement of the calculations 
Also, in many applications it will be found that the exponent m will be prac 
cally the same for all pipes in the system. 


Exponential Equations and Complex Pipe Systems 


In complex systems, pipes may occur in series, in parallel, or in net- 
works. The exponential forms of the resistance relations are especially co 
venient for solving problems of these types. 

Consider n pipes in series, each with the corresponding resistance 
relations 


cape 
h = KreaQ) 
ba ©) 
2 
hy = K, Q 
2 2 ( 
§ 
™m 
hy = K Q 3 
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4 the total head loss for the n pipes in series is 


) j=n m) m, m 

; Zh, = K\Q + K_Q oS oa Adega a a (11) 
Sue ee ° ° ° 

| j=L5j 1 2 n 


it is possible to choose exponential relations such that the actual resistance 
naracteristics of each pipe can be closely approximated using the same ex- 
pnent m; that is, if 

| 


Ee RR Re Ne ea: aan 


,en the total head loss reduces to the simple form: 


— . *. . . am 
Zh, = i ‘3 oe: rt ros oe f ee 
c t 
jen jane. 
Bihar] 2. KL Or 
j=) ej j=l j (12) 
ges. 
he quantity | z= K,,| may be thought of as an equivalent resistance 
ie ae 
ctor Ko so that the total head loss could be expressed as 
s 
jen ee 
De ee (13) 
t=L s 


se of the constant exponent m makes this analogous to the “equivalent pipe” 
ethod and subject to the same limitations referred to previously. However, 
e resistance characteristics of each pipe is characterized by the value of 
9 rather than by some fictitious length or diameter. Equation (10) permits 
e direct solution for either the total head loss or the discharge. 

If n pipes are arranged in parallel so that the head loss is the same for 
ich pipe then the total discharge through the system is the sum of the dis- 
larges through each individual pipe or 


jen 
= - ae 14 
poFet Oey Syria Fatt on (14) 


ty hy 1/m é 
it since for each pipe h,;=K,Q” or Q= = the total discharge 


°o 


r may be expressed as 


; h, \ 1/m, he 1/m, hy 1/m, (15) 
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in which hf, = hfo 5: phe he, . Again if it is possible to choose an exponentia 


resistance relation so that m is the same for each pipe, a simplification re- 
sults giving , 
1/m 1 1 1 


Q,=h + $0 «e+ 
Peet x i/m” x Tm Ro coe (1 
oi 2 n 
or 
t m 


bi item RIES m (1 
K K 
Pa 25 °n/ 


Thus an equivalent resistance factor Ko, may then be defined for the syster 
of parallel pipes as 


ae. 1/m 1/m ; 
Cia Me 
OF 22 n 


(1 
1 
FE ME PRS Sean 
j=l K, 
The head loss for the system of parallel pipes may then be expressed as 


Pp 


If the flow through any one pipe in the system is desired it may be found by 
equating the head loss in the particular pipe to the head loss through the 
parallel system 


or 


Many complex pipe systems may, for the purpose of calculations, be con 
sidered as combinations of parallel units and series units so that the entire 
system can be analyzed by repeated application of the relationships for 
parallel pipes and series pipes. Thus if the Kp values are found for each~ 
pipe in an exponential relationship with a constant exponent, equations (10), 
(11), (16), and (17) will permit a complete solution in terms of the Koy value 
for the individual pipes. 

A complex pipe system which includes some pipes which are common to 
adjacent loops can not be considered as a combination of parallel units and 
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ies units and must be analyzed as a network. Although the methods for 

lysis of networks will not be considered here, it is pertinent to note that 

alysis. resistance relations will be most convenient for use in this type 
alysis. 

; 

| Values of the Exponent m for Water Mains 

Although the Hazen-Williams equation which uses an exponent m of 1.85 

commonly used for the design of water mains in city distribution systems, 

it results appear to indicate that this is more representative of new pipe 

m of pipe that has been in service some time. It appears that the increased 

face roughness on the interior of the pipe not only increases the resistance 

ving a lower Hazen-Williams C value) but also requires a change in the 

xonent m to a value more nearly equal to 2.0. 

‘Lamont(14) presents resistance data selected from more than 200 records 

resistance measurements. In Fig. 2, the results from ninety-seven ex- 

"iments on new pipes of all types are plotted on the modified Moody dia- 

um indicating that many of the results lie in the transition range between 

ooth and rough conditions with the exponent m varying from 1.80 to 1.90. 

Fig. 3 is a similar plot for the approximately 80 records, Lamont re- 

‘ted for “old and slimy” pipes. These results clearly indicate that a larger 

ue of the exponent m is required for these pipes than for new pipes. For 

st of these records an exponent of nearly 2.0 is indicated. 

Recent tests on water mains in Austin, Texas, when plotted on the modi- 

1 Moody diagram have also indicated that values of the exponent m are 

rly 2.0. These results are summarized in Table III. 

The exponent of nearly 2.0 obtained from resistance tests on old pipes 

icates that for pipes that have been in service for some years the resist- 

-e characteristics could better be represented by an equation of the 

nning type or the Darcy-Weisbach equation with a constant f than by the 

umonly used Hazen-Williams equation. 


CONCLUSIONS 


The modified Moody diagram with lines of constant exponent m is a con- 
ient aid in determining the exponential type resistance formula applicable 
1 particular pipe for a limited range of Reynolds numbers. 


An infinite and continuous family of exponential formulas with the exponent 
ging from a minimum of about 1.70 to a maximum of 2.0 is required to 
rectly represent all pipe flow conditions. This explains why such a multi- 
e of exponential formulas have been presented in the literature on pipe 
istance. 


The commonly used exponential formulas and the equivalent pipe concept 
yased on these formulas are sufficiently accurate for many practical ap- 
ations. With some combinations of pipe sizes, roughnesses, and flow 

28, resistances may occur which cannot be explained by the common ex- 
ential formulas but which can be understood in terms of the concepts here 
sented. 
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Resistance Records on New Pipes 


4. Exponential formulas as developed from the modified Moody diagram fo 
particular system of pipes may be used conveniently in the analysis of ser: 
or parallel pipes or a combination of these. 


5. Results of resistance tests on water mains as reported herein indicate | 
for new pipes an exponent m from 1.80 to 1.90 is common which is in fair 
agreement with the Hazen-Williams formula, however, for old pipes an ex: 


ponent of nearly 2.0 is common which is in closer agreement to the Darcy 


Weisbach or Manning equations. 


A 
Cc 


List of Symbol 


Cross sectional area of flow passage - sq. ft. 
Coefficient in the Chezy equation 
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Resistance Records on Old Pipes 
) Coefficient in relation between f and Np 
i Coefficient in Hazen-Williams equation 


Inside diameter of pipe - ft. 

Inside Diameter of pipe - x inches. 

Friction factor in Darcy-Weisbach equation. 
Gravitational constant - ft/sec? 

Head loss due to pipe resistance - ft. of fluid 
Summation index ; 
Coefficient in equation hg = KQ™ 


Resistance factor a coefficient in equation he = K,Q™ chosen to be 
consistent with Moody diagram and Darcy-Weisbach equation 


ee se cee ee. 
< 
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TABLE III 


Resistance Test in Water Mains at Austin, Texas 


Pipe Age Diameter Friction Reynolds e e Export 
Description Yrs. inches Factor Number D ft. ne 
f me 
x10 


ee 


Cast Iron - 


Coal Tar Coated 18 24.0 0.0467 5.15 0017 5-0; 034 ae 
Cast Iron - 

Cement Lined 3 24.12 0.035 6.10 0.0075 0.015 Za 
Cast Iron - 

Tar Coated 45 20.00 0.036 3.96 0.008 0.013 yA. 
Cast Iron - 

Cement Lined 5 24.12 0.026 9.70 0.003 0.006 Za 
Cast Iron - 

Unlined ial 24 oe 0.041 9.20 0.093-70)026 2a 
Cast Iron - 

Coal Tar Lined 27 Wer 2 0.078 4.14 0.06 0.061 Zz 
Cast Iron - 

Cement Lined 9 6.0 Og OZ i236 0.001 0.0005 be 
Cast Iron - 

Coal Tar Lined 18 8.2) 0.082 2.29 0.065 0.045 28 
Cast Iron - 

Cement Lined 6 20rd 0.020 1.65 0.0008 0.0014 1 


Ky Equivalent resistance factor; Ko, for series pipes; K\,_ for parall 
pipes °p ‘ 


L Length of pipe - ft. 
Exponent in equation hg = Ko 
n Roughness coefficient in Manning equation; number of pipes in seri¢ 
or parallel. a 
Q Volumetric rate of flow - cfs : 


: 
R Hydraulic radius - ft. 


PIPES 


Reynolds number = ve 


Slope of piezometric head line = : 
Mean velocity = : ft/sec 


An exponent 

Linear measure of absolute roughness - ft. 
Specific weight of fluid - lbs/cu ft. 
Dynamic viscosity of fluid - lbs sec/sq ft 


Kinematic viscosity of fluid - sq ft/sec 


APPENDIX 


The following examples illustrate the steps in obtaining an exponential 
lation with the aid of the modified Moody diagram. 


& 


. 


. 


Find the exponential equation for 1000 ft. of 12-inch diameter steel 
pipe (e = 0.00015) for velocities of the order of 4.5 ft/sec with water 
at 60° F. 


Seep aaGOL Gp 
Dp URS : 
4.5 (1 piatyasiges 


R % 
L<2i 10 > 


a 
" 


From the modified Moody diagram f = 0.0155, m = 1.88. 


Mie 0.0155 (1000) (42500 ae ae 


2g (1) 2g 


nis 7 pv = Z(1) (4.5) = 3.53 cfs 


be 
K, = 2s epi Betas yeteey 
Q 3,53 1.88 
jpiwet agit, wise 


Find the exponential equation for 100 ft. of 6-inch cast iron pipe (e = 


_ 0.0009) for velocities of the order of 5 ft/sec with water at 60° F. 


ees =r) OROOL8 
D 


39 


40 


3. 


oi 


March, 1959 


Ne To ace ee 2.06 x 10° 


1.21 %.30°" 


f=. 00235; re =. oo 


2 
pow 
af we e182 ft 
b= £ & 3 
Qa=f (0.5)* (5) = 0.98 cfs 
h 
K, = = tae = 188 
ao wine) ee 
bh, =p1.88Q i?” 
£ 
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CHARGE CHARACTERISTICS OF RECTANGULAR THIN-PLATE WEIRS@ 


: 

| Closure by Carl E. Kindsvater and Rolland W. Carter 

: 

SARL E. KINDSVATER, 1 M. ASCE and ROLLAND W. CARTER,?2 A. M. 

‘= —In concluding this discussion of a paper on the subject of weirs, it 

ms appropriate to observe that few subjects in the literature of hydraulics 
more revealing of the personalities of its students. It is indeed remarka- 
that a comparatively simple physical phenomenon has inspired such an 
unding number of printed words, formulas, laboratory investigations, and 
wgreements. Over a period of 250 years, nearly every name in the history 
ydraulics is identified with contributions to the great store of weir litera- 
>. It is not amazing, therefore, that the most common ingredient in this 
sritance is repetition. 

Much more discouraging than the repetition which characterizes weir 
rature is the fact that some of the best thought on the subject is buried in 
avalanche of printed words. On the other hand, the substantial area of 
rance regarding weirs which still exists is often dismissed with the casu- 
mplication that the solution to the problem is contained in a given formula 
hich is then quoted without reference to the qualifications which were pre- 
ibed by its author. 

ride and prejudice seem to have played a large part in the drama of the 
ple weir. Counter-claims of accuracy and scientific truth have been 

ed across battle lines drawn up behind purely empirical formulas based 
imited, often questionable experimental data. The futility of one such 
iment, which already occupies many pages in the Transactions, is revealed 
‘ig. 4(a). Here it is clearly demonstrated that, regardless of the quality 
azin’s data, the form of discharge equation used by Bazin and others is 
lamentally inadequate. 

teturning to the present, the writers are grateful for the constructive dis- 
sions of their paper which were submitted by Messrs. Sarpkaya, Kolupaila, 
ell, Paull, Oki, and Carstens. It is regretted that several correspondents, 
iding some respected European hydraulicians, were unable to submit 

r discussions within the time allotted for discussion. 

‘he writers share with Mr. Sarpkaya his concern with the practical limi- 
ms of empirical equations. They are aware of the dangers of extrapolating 
irical relationships beyond the range of conditions represented by the 


ro . Paper 1453, December, 1957, by Carl E. Kindsvater and Rolland W. 
Beier. 

tegents Prof., School of Civ. Eng. Georgia Inst. of Technology, Atlanta, 
ja., and Consultant, Geological Survey, U. S. Dept. of Interior, Washington, 
aC. 

thf., Research Section, Surface Water Branch, Water Resources Div., 
‘eological Survey, U. S. Dept. of Interior, Washington, D. C. 
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experiments from which they were derived. Persons who publish empiricé 
equations are obliged to state very clearly the purpose and limitations of 
their work. It is the reader’s obligation to observe these limitations in the 
application of the results. 

The writers have made it quite clear that they were concerned with the 
weir as a discharge-measuring device, and that, for practical reasons, the 
verification of the proposed formula was limited to “water at ordinary 
temperatures”. Furthermore, the range of values of each variable involve 
in the experimental verification is explicit in the original paper. Without 
apologizing for the limitations which time and money placed on their work 
the writers are, like Mr. Sarpkaya, hopeful that others will extend its bour 
ries. They are inclined to observe that, apart from the acknowledged limi 
tations in their study of the fluid-property variables, their work covers a 
broader range of geometric conditions than any other single work known. 

It should be quite clear from the original paper that values of k}, and ky 
are dimensional quantities derived from experiments with a single liquid : 
a practical range of geometric conditions. The general relationships betw 
kp, kp, R and W are carefully detailed. One accustomed to using the Reyn 
and Weber numbers to describe the relative influence of viscosity and sur 
tension, respectively, would not expect “glucose, honey and molasses” to 
exhibit unique characteristics as long as they behave as true fluids. 

Mr. Sarpkaya’s Eq. (3) is a general statement of the slope-intercept fo: 
of the Ce curves (writers’ Fig. 10), which are believed to be descriptive o 
the interrelationship of the principal geometric variables. His Eq. (2), on 
the other hand, involves a dimensionless representation of the fluid-prope 
adjustment factors, ky and kp. It is not clear why he presumes that an ad 
ditional fluid-property adjustment should be made to the relationship betw 
the geometric ratios. Indeed, Eq. (5) appears to hold little promise of sin 
fication or clarification. And, while it is considerably more complex than 
writers’ equation of discharge (Eq. (11), it is certainly no less empirical. 
Subsequent algebraic manipulations and the unlikely assumption that kp is 
independent of surface tension do not lend credence to Mr. Sarpkaya’s pro 
posed equation. 

Mr. Kolupaila incorrectly attributes certain objectives to the writers ( 
“a search for higher precision [than + 2% ] was the purpose of the new inv 
gation”) and their sponsors (e.g., “The sponsor . . . expecteda. .: formu 
correct to + 5%”) He is also incorrect in attributing to the writers’ thoug 
which are not to be found expressed in the original paper (e.g., “Their opi 
is based on usually assumed values of & of about 1”, and “the conclusion ¢ 
the Walchensee tests was expressed as being unfavorable by the authors”) 

The writers agree with Mr. Kolupaila that the flow pattern for rectang\ 
weirs might be subject to complete mathematical analysis if “a certain la 
streamlines” were known. Unfortunately, the solution is hardly forthcomi 
in view of his subsequent qualification that such a law remains to be “es- 
tablished or approximation assumed” along with other “small correction 
factors”. Mr. Kolupaila’s remarks concerning kj and & are somewhat co 
dictory, but, in general, they appear to substantiate the writers’ conclusi 
regarding the tenuous relationship between these two coefficients. a 

Mr. Kolupaila makes reference to many foreign publications on weirs, 
most of which had been read and evaluated by the writers. Failure to ref 
certain of these publications in the original paper is evidence that they we 
not considered of particular value. It is not clear why Mr. Kolupaila refe 
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e foreign publications without identifying and supporting his own prefer- 
c€é among the diverse opinions and formulas proposed therein. A casual 
view of recent European reference works indicates that few of the papers 
erred to have been recognized as substantial contributions in their own 
elands. Finally, an examination of the paper which describes the 
lenense tests will show that Mr. Kolupaila, not the writers, must have 
sinterpreted the conclusions drawn by Kirschmer and Esterer. The follow- 
¢ free translation of their conclusions clearly supports the conclusions 
= by the writers: 


“The results of the weir measurements substantiate the criticism of 
|this method which has been expressed by others. Indeed, general 
‘reliance upon one of the rational formulae has not yet been satisfactori- 
ly proven”. 


‘Mr. Powell emphasized the unreliability of coefficients derived from ex- 
riments on physically small values of P. The writers have acknowledged 

.§S uncertainty as well as the uncertainty of the results of tests made with 
ry small values of b and h in their investigation. 

Mr. Powell’s suggested method of correcting for the distribution of veloci- 
in the approach channel is a very worthwhile and practical contribution. 

ie results which he obtained from an analysis of the Schoder and Turner 

ta are encouraging. It is hoped that additional tests made to define the “r” 
ction will refine the method. Thus, the weir which is unavoidably located 
a channel which produces a non-uniform velocity distribution might still 
Salvaged as a satisfactory meter. As pointed out by Mr. Powell, the ratio 
- Va/Vp is essentially a characteristic of the installation which, like the 
nensions of the weir, need be determined only once. It follows that his 
thod of adjusting the coefficient of discharge is completely practical. It 
ems reasonable to assume that r is a function of h alone. Mr. Powell’s 
lingness to settle for “errors of the order of one percent” is quite 

alistic, if not optimistic. 

The writers appreciate Mr. Paull’s kind remarks. They regret that they 
inot reciprocate with an offering of dependable information regarding the 
ccuracies to be expected under a full variety of applications. Disregarding 
ch very real problems as variations in velocity distribution and “non- 
ndard” geometric conditions, the order of magnitude of the errors to be 
ected are suggested by the scatter of the plotted points on the graphs show- 
‘the writers’ data. — 

Mr. Oki brings to the discussion some new experimental data and formulas 
viously unknown to the writers. They are pleased to acknowledge his 
erest and criticism. In his Fig. 1 Mr. Oki proposes to demonstrate a re- 
ionship between C and h/b. In their original paper the writers state, “The 
ependent influence of this ratio (b/h) is believed to be negligible over the 

| practical range of the other variables. An earlier investigation at the 
orgia Institute of Technology supports this conclusion. . . . A few recorded 
rts to incorporate the b/h ratio in discharge formulas are believed to be 
ed on misinterpretations of influences related separately to the magnitudes 
) and h.” 

The writers believe that this statement is actually supported by Mr. Oki’s 
. 1. Indeed, closer study will reveal that the slope of the curves on the 

ire shows the effect of b, whereas the vertical distance between the curves 
ws the effect of h. In fact, all of the data on Fig. 1 can be replotted to 
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show a relationship between Cg (= C + 2/3 V2g) and h for various constant 
values of b. However, if this is done, Mr. Oki’s data appear to show a con- 
siderably larger influence due to both h and b than has been reported by the 
writers and others. 

The writers are especially puzzled by some of the very low values of C, 
shown on Mr. Oki’s Fig. 1. All values of Cg for h/b greater than 0.5, for e 
ample, are less than the “theoretical” minimum, 0.611. This trend is be- 
lieved to be unique in the literature. It is implied, of course, that the de- 
crease in Cg is related to the increase in h/b. The writers can only repeat 
the conclusion, based on their own work, that the h/b ratio is not believed t 
be significant. In support of this conclusion they submit the observation thi 
the data presented in the original paper cover a range of values of h/b fron 
0.03 to 6.0. 

Mr. Carstens lists several reasons for the difference in the results ob- 
tained by applying the writers’ method of analysis to the experimental data 
obtained by, different investigators. The writers agree with his general dia 
nosis, but they fail to see the logic in his conclusion that the maximum exp 
mental error in one set of data is an explanation of the difference between 
curves which were drawn to represent mean values for different sets of da 
Such an unusual conclusion ignores, for example, the systematic variation 
with velocity distribution, weir crest condition, etc., shown by Schoder and 
Turner and others. To make their own position clear, the writers go on 
record as pointing to Mr. Carstens’ item (b), namely, the “geometry of the 
experiments”, as the major cause of the difference between the results of 
different experiments. A large part of the tremendous store of weir litera 
ture is evidence to support this conclusion. 

Mr. Carstens’ light dismissal of the “small differences” between variou 
formulas ignores the purpose of most of the work done in this field since tl 
original derivation of Eq. (15) by Weisbach. His comparison of the “theore 
cal” and measured values of Cg is a classic demonstration of circumstanti 
evidence rather than theoretical proof of the merit of the orifice analogy. 
The writers are pleased to endorse his conclusion that too many people wr 
empirical equations (or equations that are too empirical !), but they fail to 
the relationship between this popular thesis and the statement quoted from 
the original paper which preceded his pronouncement. They are also move 
to observe that the weir which is the subject of Mr. Carstens’ attention (i.¢ 
b/B = 1.0) is only a limiting case of the problem under investigation, that 
values of h/P greater than 2,0 are impractical for reasonably accurate flo’ 
measurement, and that the fluid-property effects associated with small val 
of h and b cannot be ignored for all practical purposes. 

In general, the discussers did not challenge the writers’ candid disresp 
for the traditional “theoretical” equation and the orifice analogy. However 
most of them appeared to favor the use of a coefficient Cg which is equal t 
the writers’ C divided by a quantity (2/3 2g ) derived from the orifice 
analogy. The principal advantage attributed to Cq is that it is dimensionle 

The writers are staunch proponents of dimensionless ratios. A review 
their publications will substantiate this statement. As engineers as well a 
researchers, however, they strive to be practical. In the present instance 
they question the logic of straining to make the coefficient of discharge di- 
mensionless when the equation of discharge contains other terms which ar 
dimensional. In the writers’ equation, for example, evaluation of the disck 
depends on knowledge of ky and kp, quantities which are absolute measure: 


'3 DISCUSSION 49 


» combined influences of viscosity and surface tension. Similarly, in Eqs. 
) (Bazin), (19), (Rehbock), and (20) (S.1.A.), dimensional terms which are 
ated to the magnitude of h restrict the application of those equations to 
ter at ordinary temperatures. 
‘It has been observed that the quantity 2/3 y2g is obtained from the deri- 
ion based on the orifice analogy. Bazin and others have at one time pro- 
se equations in which 2/3 y 2g was replaced by ye or Ve. In these alter- 
e forms of the discharge equation the coefficient is quite dimensionless, 
t the writers are loath to use them because of the almost exclusive division 
preference between C (widely preferred by engineers in practice) and Cq 
referred by some researchers) as defined herein. While they find it diffi- 
lt to get excited about the choice, the writers decided to use C mainly be- 
use it concedes to the practicing (American) engineer the margin of con- 
nience. 
The writers freely admit the advantage of C, (rather than Cg. See Eq. 
5)) for certain purposes. A reason seldom mentioned, for example, is that 
. is remarkably similar for different shapes of notch weirs. Thus, the re- 
archer finds it convenient to use Cc; to correlate meagre data from differ- 
t sources. However, this reason does not justify the preservation of the C, 
rm of the discharge equation for practical use. It is hoped that this brief 
scussion of the relatively insignificant difference of opinion regarding C, 
j and Ce does not detract from the purpose of the original paper. 


Corrections.—On page 5 of Proc. Paper 1453, in the third line following 
. (6), the term V¥o/ p should be changed to Yo/p. On page 27, in the first 
.e of the paragraph containing Eq. (22), the word “new” should be changed 
“few”. F 

On page 1690-24, line 5, there should be no parenthesis after the 3.23. 
On page 1690-25, eighth line below Eq. (F), “which Eq. (D)” should be 
yhile Eq. (D)”. 

The term V2z in the third line of the fourth paragraph on page 1690-22 
ould be 2/3 Ve. 


ae a ee Pa 
= 


may WI 9cl # 
¢ © 


oO x ous% 49 
*- = 


51 


THE TOTAL SEDIMENT LOAD OF STREAMS2 
Closure by Emmett M. Laursen 


EMMETT M. LAURSEN, | A. M. ASCE.—At first reading the three ways of 
ying the problem of total sediment load enumerated by Messrs. Garde and 
ertson seem wise, considered and reasonable. Repeated reading, however, 
is one to doubt, and finally reject, most of what is said or implied in this 
tion of their discussion. Not enough is known about fluid flow to theoreti- 
y describe in detail the state of affairs at a complexly rough fixed bounda- 
How then does one even begin a “theoretical” evaluation of the total sedi- 
wt transportation rate? The total load is undoubtedly the sum of the sus- 
ded and the bed load but which bed load equation does one choose and how 
S one evaluate the reference concentration of the suspended load except by 
itrary assumption? (For a detailed discussion of the Einstein procedure 
the modified Einstein procedure the reader may be interested in Refer- 
e (6)). This leaves dimensional analysis and intuition as possible methods. 
1ensional analysis can only lead to the grouping of chosen variables. Ergo, 
ition must reign supreme. Unfortunately, this state of affairs is all too 
», Yet intuition is not entirely free, but must be bound by what is known. 
\s presented in the discussion, the relationship for the total sediment load 
forward by Messrs. Garde and Albertson seems to have been developed 
tly by intuition using dimensional analysis to form the dimensionless pa- 
eters and plotting and curve fitting to obtain the relationship between the 
ameters. In the form presented it appears to be different from other sedi- 
it transport equations. However, its similarity to most other equations 
ymes obvious if it is rewritten. Using the symbols of the original paper 
Yo; V, = V70/P /2 
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» that V is a function of d alone, and the curve on the log-log plot can be 
-oximated by the straight line Y = 0.12 d-3/2, This equation can then be 
her transformed through the use of the relation qg = const. ct q and the 
ning equation to read 


3 yt 


=C n ah os 
eae qt Yo 
re the coefficient C has been allowed to absorb the various factors which 


constant or almost so. Similarly transformed, and with the critical term 
ped, other equations are very similar.\1 


roc. Paper 1530, February, 1958, by Emmett M. Laursen. 
.ssociate Prof., Dept. of Civ. Eng., Michigan State Univ., East Lansing, 
Nich, Zang 
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Note that in some of these equations a part of the effect of the sediment sii 
is contained in the coefficient. : 

Mr. Bogardi’s interesting contribution, like that of Messrs. Garde and 
Albertson, is more a reporting of similar work than a discussion of the co 
cepts of the original paper. The relationships presented delineating the be 
configuration could be very useful if confirmed over a wider range of con- 
ditions. The writer hopes that in Mr. Bogardi’s projected paper, the ratio 
of at least the parameters, and preferably also the relations, will be dis- 
cussed at greater length. The use of the total tractive force (or the shear 
rived therefrom) in the bed stability factor and the diameter of the seceding 
by itself seem to be questionable procedures, 

The discussion of Mr. Bondurant is of the type that the writer had most 
hoped for—the application of the proposed relationships to field measure— 
ments. Considering the change in scale and general environment between 
laboratory flumes, which provided the data upon which the relationships w 
built, and the Missouri and Arkansas rivers, the factor of four, or the shit 
in the {(U,/w) curve, is not too disheartening. The discrepancy could be ¢ 
either to this gross change in scale or to deficiencies in the proposed re- 
lationships, or possibly to both reasons. The fact that the original curve ¢ 
scribed the Niobrara measurements so well and that the shift of the curve 
the Missouri river data from Omaha was less than for the data at Kansas. 
would indicate that some factor such as the Reynolds number might be usée 
in collapsing the different results. More field measurements might perm: 
the development of such a stop-gap modification and thereby enhance the - 
usefulness of the proposed relationships. 


TER DISTRIBUTION DESIGN AND THE McILROY NETWORK ANALYZER@ 


| 
| 
| Closure by M. B. McPherson and J. V. Radziul 


M. B. McPHERSON, ! A. M. ASCE and J. V. RADZIUL.2—The supple- 
ntary experiences, observations and comments contributed by Messrs. 
wes and Branscome, Mr. Cole, and Mr. Lomax,? are sincerely appreciat- 


The need expressed in the paper for publication of a proven standard set 
nstructions for water distribution network analysis for a stored-program 
2 of computer has been satisfied by Messrs. Graves and Branscome. 1-a) 
‘ir paper and the authors’ appeared in the Proceedings simultaneously. 
sequently, Mr. Hamblen(1-c) offered a program suited to calculations for 
yer networks. With the publication of these programs there is no longer 
alid excuse for rough approximation by means of extensive shortcuts in 
design or investigation of complex water distribution networks by consult- 
engineers and utilities simply because a McIlroy Analyzer does not happen 
e readily available to them. 

Mr. Lomax says, “The possibilities of the digital computer and the ana- 

sr have not been fully developed,” and illustrates some of these potentials 
temizing several unusual problems he has analyzed. The authors are in 
agreement with his statement: “The combination of the analyzer and the 
yputer (at his institution) will provide complementary rather than competi- 
Services,” but also foresee valuable collateral uses for both devices in 
ribution system analysis by the larger utilities. 

ince large networks are almost always designed on the basis of peak day 
jitions, “computations on the basis of the largest rates of demand” 

aves and Branscome) are normally needed regardless of the device used 
erforming the computations. Without equalizing storage, the demand of 
peak hour of the maximum day normally constitutes the critical design 
lition for large networks. With storage, the minimum hour of the maxi- 

n day demand is apt to be more critical because the storage is necessarily 
g replenished at a high rate at that time and because the flow to storage 
nally traverses a long distance across the network. The design of arterial 
ng improvements is a trial and error procedure. When estimated future 
ands (always taken at rates greater than current demands) are being 

ied, the authors use an initial arterial network which is usually a moder- 
modification of existing piping but which will probably yield excessive 
llosses. They have been able to reach a more satisfactory final piping 


Sroc. Paper 1588, April, 1958, by M. B. McPherson and J. V. Radziul. 
tesearch Engr., Philadelphia Water Dept., Philadelphia, Pa. 

Vater Distribution Design Engr., Philadelphia Water Department, 
hiladelphia, Pa. 

Miscussions appeared in Proc. Paper 1856, p. 85-88, Nov., 1958. 
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arrangement more quickly in this manner than by starting with a network : 
which would undoubtedly prove to be more than adequate. Hence the stater 
in the paper: “--- needed system revisions are evident at inflow rates wel 
below the design rates; obviously needed piping improvements are then mai 
prior to any formal or complete test runs.” 

Messrs. Graves and Branscome correctly note, with reference to Fig. 6 
of the paper, that the statement: “A revised storage rate requires a new 
storage curve,” is not clear. The balancing of demands, storage rates and 
near-uniform pumping rates is difficult to describe without illustration. Hi 
ever, a detailed example of balancing procedures using hour-by-hour analy 
for the City of Philadelphia’s Belmont High Service District, with a single 
pumping station and one storage site, is to be published soon. 2 

Mr. Cole raises the question of economic justification of the equalizing 
storage required for near-uniform pumping versus the lesser storage need 
for stepped pumping. Comparisons are extremely complicated and no uni- 
versal rule can be applied even between service districts of the same wate 
works. The principal benefit sought in the installation of equalizing storag 
is the provision of more stable pressures and less service interruption tha 
can be achieved by direct pumping. Unfortunately, the storage volumes ne¢ 
ed to secure near-uniform pumping, augmentation of pumping station outpu 
during fires, and outage relief, are neither a summation of the volumes net 
ed for each purpose nor a readily determined overlapping of volumes. 

Consider first, peak day storage. As an example, if the maximum day 
demand is about 150% of the annual average, 10% of the average day demai 
is needed for near-uniform pumping on the average day and 10% of the mai 
mum day demand is needed for near-uniform pumping on the maximum day 
the storage volume required for the maximum day is thus about 1.5 times} 
volume required for the average day. This theoretical ratio ranges betwee 
1.3 and 2.1 for the five districts cited in Fig. 4. The ratio is undoubtedly 
higher in many systems across the country. Maximum day storage raises 
and stabilizes critical low pressures under non-emergency operating con- 
ditions. Benefits received on days of lesser demand might be regarded as 
bonuses not chargeable to the cost of the maximum day storage. If the sto: 
provided is not sufficient to meet peak day variations, the resultant higher 
pumping rates on peak days will generally require a more expensive distri 
bution network than would be satisfactory for near-uniform pumping with n 
storage. ; 

At what period during the design maximum day should fire requirement 
be met? Among other things, the N.B.F.U. rates a district on ability to de 
liver the required fire flow during a district demand equal to the maximun 
day rate. Philadelphia’s design policy is to provide a working elevated 
storage volume (upper 25-feet) adequate to meet future maximum day dem 
under hypothetically perfect uniform pumping (exceeds somewhat the quant 
needed for near-uniform pumping). An additional, separate increment of s 
for fire flow augmentation is not provided. The fire flow increment is to k 
delivered from the pumping station. Special design analyses are made to 
insure adequate fire pressures under these conditions. In large districts, 
flow and pressure requirements for fires seldom dictate additions over an 
above the piping needed to satisfactorily meet peak day demands. 3 

For what duration and for what hour of the day should outage relief b k 
vided? In Philadelphia it is not considered economically feasible to earm: 
a special storage increment, or special piping, for outage relief alone. _ 


> 
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ever, the extra benefit of added outage relief has been a deciding factor 
iepting some otherwise marginally justified, needed improvements. 
niess the total equalizing storage volume provided is at least sufficient 
eet average day demand variations, the benefits which will be derived 
the storage on a peak day will be almost inconsequential. Therefore, 
inimum amount of storage worthy of serious consideration should be the 
me required to meet average day demand variations. With this mini- 
m, near-uniform pumping may often be feasible for more than half of the 
-’s demands, as inferred in Fig. 3. The annual cost of pumping will be 
atly influenced by the large share of demands thus met by near-uniform 
ping with a single pump. The cost of an increment of storage above the 
imum to meet peak day fluctuations can be readily compared against the 
- of additional piping which would be required to provide similar service 
ssures by means of direct peak day pumping without storage. The more 
ortant comparison, from a cost standpoint, between the use of minimum 
‘age with stepped pumping on the peak day and either direct or near- 
orm pumping would be purely speculative. For stepped pumping, it is al- 
st impossible to realistically predict the most probable pattern of peak 
demands and the sequence and timing most likely to be followed by the 
ion operator in meeting these variations. 
The dollar value of the more reliable overall service afforded by equaliz- 
storage is largely intangible despite its obvious importance. Unfortunate- 
he cost of an increment of storage especially set aside to more adequately 
t fire demands cannot be equated in direct terms to a specific expected 
iction in fire insurance policy rates. 
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SYNTHETIC FLOOD FREQUENCY? 
Discussions by Ven Te Chow and J. L. H. Paulhus 


"EN TE CHOW,! A. M. ASCE.—The author has presented an interesting 
redure for synthesizing flood frequencies through the use of a rational 
qula and unit hydrograph. In essence, the lengthy derivation of the pro- 
ire may be boiled down to the effect that the flood discharge can be com- 
d by the conventional rational formula and then reduced in the magnitude 
rder to convert the value for the rainfall frequency to that for the runoff 
1¢ Same frequency. When the conventional rational formula is used, it is 
cal to accept that the frequency of the computed discharge is compatible 
the frequency of the rainfall intensity which is used in the formula. Ow- 
to the detention and retention characteristics of a watershed, the actual 
uency of the resultant discharge is in general much higher than the 
uency of the corresponding rainfall. If the computed discharge takes the 
é frequency of the rainfall, the computed value should be therefore re- 
od. 

‘he above discussion may be expressed by a modified rational formula: 


Q = KCIA (14) 


re Q is the peak discharge in cfs., K is a reduction factor, C is a runoff 
ficient, I is the rainfall intensity in in. per hr., and A is the drainage area 
ores. 

‘he use of the unit-hydrograph theory in the author’s derivation indicates 
the rational formula in the form of Eq. (10) is a special case of the unit- 
Ograph method, A elie mathematical proof of this point has been 
onstrated by Nash. (1) 

q. (14) is in fact given in the form of Eq. (13) in the paper. The author 

d that the value of K in both formulas may be considered as a constant, 
is equal to 0.44 for overland flow on turf and equal to 0.77 for natural 

S, sewered areas, and overland flow on pavements. In other words, the 
tionship between the rainfall and runoff for a given frequency is assumed 
tant. It is believed that the author has accepted that rains of the same’ 
tion and frequency occur with associated runoff producing conditions that 
from time to time depending on antecedent factors. It may be then 

oned that an averaged rainfall-runoff relationship for that group of rain- 
would produce an acceptable runoff volume-frequency value. The 
fulness of this reasoning can only be verified by the use of actual data. 
ugh the author has shown several actual examples for this purpose, it 


roc. Paper 1808, October, 1958, by Franklin F. Snyder. 
rof. of Hydr. Eng., Univ. of Dlinois, Urbana, 11. 
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would be interesting to see how this reasoning would apply to many other 
areas which are different from those investigated by the author. 

It is generally recognized that two major factors that lead to runoff are 
the storm rainfall, a climatic factor, and the antecedent conditions, a physi 
graphic factor. The probability occurrences of rainfall and runoff are well 
known, while the probability occurrence of the physiographic factor is a 
subject that has not yet been explored. From a probability viewpoint, how- 
ever, the physiographic factor is a statistical variable just as the rainfall ; 
runoff are statistic variables.(2,3) The frequency of runoff is the product ¢ 
the frequencies of rainfall and physiographic condition. When an average 1 
lationship between rainfall and runoff frequencies is assumed, the frequenc 
of physiographic variation is automatically taken as constant. Whether thi: 
is true or not remains to be seen from future verification by data of variot 
geographical characteristics. 
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J. L. H. PAULHUS, ! M. ASCE.—A novel and generally logical approach 
the flood-frequency problem has been developed by Mr. Snyder. His pro- 
cedure appears suitable for deriving flood-frequency data for very small 
areas, especially if they have a large percentage of impervious surface, p: 
vided representative rainfall-frequency data are available. There are, hor 
ever, several weak points in its development that make Snyder’s approach 
less adaptable to larger basins. In general, when representative rainfall- 
frequency data are used, the weaknesses would indicate a tendency for the 
procedure to yield results somewhat higher than might reasonably be ex- 
pected. 

Perhaps the biggest drawback to the successful application of Snyder’s 
cedure, even when properly developed, is the lack of means for selecting t 
most representative rainfall-frequency data to use. Snyder advises that v. 
ation in annual precipitation must be considered in the selection or adjusti 
of the point rainfall record to be used but omits details. This practice doe 
not necessarily produce reliable results. Table 1 shows the different floo 
frequency values that would be obtained for Wills Creek near Cumberland, 
Md., by applying Snyder’s procedure to rainfall-frequency data, unadjustec 
except for partial-duration series, for the five nearest surrounding statio 
for which that type of data are already available. Pittsburgh and Elkins a1 
each about 75 miles from Cumberland and the other three stations are bet 
105 and 115 miles away. The two nearest stations yield widely different s 
of results, with flood-frequency values based on Elkins rainfall being abou 


1, Staff Hydrologist, U. S. Weather Bureau, Washington, D. C. 
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ste 1. FLOOD-FREQUENCY DATA (100CFS) FOR WILLS CREEK BASED ON 5 STATIONS 
| Return period (yr) 2 5 10 25 50 100 
| Obs. record, 27 yr 92 155 203 268 316 362 
| Pittsburgh, Penn. 57 78 96 119 138 159 
Elkins, W. Va. 75 WW 144 191 233 279 
. Harrisburg, Penn. 80 5 148 198 237 279 
| Baltimore, Md. 105 154 194 256 311 363 
| Washington, D. C. 125 178 225 290 358 42k 


Der cent higher than those based on Pittsburgh. Which station is repre- 
ative for the basin? Noting also that Baltimore yields results over 25 
cent higher than those of Harrisburg, with normal annual precipitation of 
59 and 36.01 in., respectively, which of those two stations, or what kind of 
istment, would one use in developing a synthetic flood-frequency curve for 
asin midway between the two cities ? 
The results of Table 1 would be changed very little by adjustment of basic 
afall-frequency data for differences in mean annual precipitation. Stations 
ind near Wills Creek Basin indicate a mean annual precipitation between 
‘ssburgh’s and Elkins’. The differences would not provide a suitable basis 
an adjustment yielding flood-frequency values approximating the observed. 
$ possible that neither Elkins nor Pittsburgh yields results comparable to 
observed because Snyder’s rainfall-runoff relation is not representative 
Wills Creek. There is a very small chance that it could be considering 
derivation, which is discussed later. 
Ine of the more serious shortcomings of the procedure is the assumption 
_ the rainfall intensity for a given return period will produce a peak- 
charge value having the same return period. For example, the 10-year 
fall intensity is-assumed to produce the 10-year flood. In this respect, 
der’s procedure is no better than formulas such as the so-called rational, 
SiA, formula. 
Annual maximum peak discharges often occur in connection with rainfall 
nsities less than the annual maximum intensities. The reason is that in 
ly regions the highest intensities occur in the seasons when soil conditions 
not generally favorable to high runoff rates. For example, Snyder would 
sumably use rainfall intensities for Washington, D. C., to develop flood- 
juency data for the Rappahannock River near Fredericksburg, Va., which 
a drainage area of 1599 square miles and a computed concentration time 
) of 24.8 hours if it is assumed that the impervious area of the basin 
unts to less than 1 per cent. The rainfall-frequency data he would use 
e obtained by the Gumbel method of analysis utilizing the annual maximum 
les for the 57-year period 1896-97, 1899-1953. The results were then ad- 
ed by Snyder for a partial duration series and converted from intensities 
mounts to yield his Fig. 2, Appendix 1. Comparison of the monthly distri- 
on of the annual floods of the Rappahannock for the total 45-year period 
3-57) of peak discharge record with that of the annual maximum 24-hour 
fall intensities used in the Gumbel analysis shows a relatively poor re- 
ship (Table 2). 
t is readily seen that the 3-month period yielding the greatest percentage 
nnual floods is March-May with 37.9 per cent. Only 15.8 per cent of the 
al maximum 24-hour rainfalls occurred in the same period. On the other 
|, while the 3-month period, July-September, yielded 57.8 per cent of the 
ial maximum 24-hour rainfalls, only 26.6 per cent of the annual floods 
irred during this season. en 
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Table 2. MONTHLY DISTRIBUTION (PERCENT) OF ANNUAL MAXIMUM FLOODS AND 
24-HOUR RAINFALLS 


Month Ann. Floods Ann. Max. 24-hr 
(Rappahannock) Rainfalls 


(Washington) 


Jan 6.7 0.0 
Feb. 2.2 Ki7 
Mar 8.9 3.5 
Apr. 17.9 5.3 
May 11.1 7.0 
Jun, 4k 8.8 
Jul. 4k 17.5 
Aug. 11.1 15.8 
Sep. 11.1 24.5 
Oct. 8.9 5.3 
Nov. 8.9 5.3 
Dec. 44 5.3 


Since the monthly distribution might reasonably be expected to be bette 
than the above for a small basin for which point rainfall would be more re 
sentative, it appeared advisable to make atest. Unfortunately, there are! 
streamflow records of adequate length for any small basin in the area of 
interest. Consequently, the writer synthesized a record of annual floods f 
the Little Falls Branch near Bethesda, Md., rey area 4.1 square milt 
using the procedure outlined in Proc. Paper 1451.(1 Hourly rainfall data 
Baltimore, Md., were used in synthesizing the annual flood record as they 
were already available to the writer. 

Bearing in mind that the writer was, at this point, concerned only with 
comparison of the monthly distribution of annual floods and annual maxim 
rainfall values, it should be noted that the accuracy of the 1-hour unit hyd: 
graph used in synthesizing the annual flood record had little bearing on thi 
monthly distribution. The reason for this is that the annual maximum 1-h 
computed runoff was almost always responsible for the annual peak disch: 
for the basin, which, according to the author, has a computed T, of 1.65 
hours. Consequently, the comparison could have been made directly by us 
the monthly distribution of annual maximum hourly runoff with no appreci 
change in the results. 

Annual maximum peak discharges were synthesized for the 61-year pe 
1894-1954, and their monthly distribution noted. The monthly distribution 
annual maximum 2-hour rainfalls at Baltimore, Md., for the 49-year peri 
1903-51, used in the Gumbel analysis for U. S. Weather Bureau Tech. Pax 
25,(2) was then determined. It was noted that annual maximum 1-hour rai 
falls had almost exactly the same monthly distribution as the 2-hour. Th 
comparison of monthly distributions of annual maximum floods and 2-hou: 
rainfalls is given in Table 3. As was expected, the latter comparison sho 
better agreement than the first. The 3-month period, J uly-September, ha 
the greatest percentage of annual floods and annual maximum 2-hour rain 
although the proportions differ considerably, 52.5 vs. 71.5 per cent, re- — 
spectively. Actually, of the 49 values of observed annual maximum 2-hou 
rainfalls, 30, or 61 per cent, produced annual maximum synthesized peak 
charges. The other 19 synthesized annual maximum peak discharges res 


from lesser rainfalls occurring with more favorable flood-producing soil 
conditions. 
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je 3. MONTHLY DISTRIBUTION (PERCENT) OF ANNUAL MAXIMUM FLOODS AND 2-HOUR 
: 
| 


RAINFALLS 
Month Ann. Floods Ann. Max, 2-hour 
(Little Falls Br.) Rainfalls 
(Baltimore, Md.) 
: Jan 1.6 0.0 
| Feb. 1.6 0.0 
| Mar. 3.3 0.0 
Apr. 6.6 2.0 
| May 11.5 6.1 
' Jun, 8.2 12.3 
i Jul. 18.0 28.6 
i Aug. 23.0 28.6 
| Sep. 11.5 14.3 
Oct. 4.9 | 
Nov. 4.9 2.0 
Dec. 4g 2.0 


“he above comparisons suggest that, other things being equal, Snyder’s 
cedure should tend to yield flood-frequency values that are too high. Un- 
unately, only one of the basins the author used has a streamflow record 
ufficient length to yield a reliable test. The Rappahannock River near 
dericksburg, Va., has a 45-year record (1913-57) of annual peak dis- 
rges. The annual floods were analyzed by the Gumbel method and the re- 
S adjusted for a partial duration series. The comparison between these 
sted values and those to be obtained by Snyder’s method are presented in 
le 4, 


le 4. FLOOD-FREQUENCY DATA, RAPPAHANNOCK RIVER NEAR FREDERICKSBURG, VA., 


(100 CFS). 
Return period (yrs) 2 5 10 25 50 100 
Obs. 45-yr record 380 590 748 960 1120 1279 
Snyder's method 384 560 720 968 1160 1360 


1 view of the author’s expressed belief that a flood-synthesizing procedure 
1 as described in Proc. Paper 1451(1) would yield more uncertain results 
his method, comparisons were made for two other basins used by him 
for which annual floods had already been synthesized. These two were 

se Creek near Leesburg, Va., drainage area 338 sq. mi., and Little Falls 
ich near Bethesda, Md., drainage area 4.1 sq. mi. The observed and 
hesized annual floods were analyzed by the Gumbel method and then ad- 

sd for a partial duration series. The results are listed in Table 5 for 
parison with the values obtained by Snyder’s method: 


le 5. FLOOD-FREQUENCY DATA, GOOSE CREEK NEAR LEESBURG, VA., (100CFS). 


_ Return period (yr) 2 5 10 25 50 100 
Obs. 26-yr. record 126 233 308 410 485 562 
Snyder's method 115 175 216 287 348 407 


Synthesized floods,55 yrs. 98 163 210 277 327 376 
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LITTLE FALLS BRANCH NEAR BETHESDA, MD., (10 CFS). 


Return period (yr) 2 5 10 25 50 100 
Obs. 13-yr. record 116 160 19] 235 270 301 
Snyder's method (Wash. data) 118 144 167 208 241 276 
Snyder's method (Baltimore data) 90 122 152 189 222 255 
Synthesized floods, 61 yrs. 100 133 158 191 217 24] 


The comparative tabulations of flood-frequency data in Tables 4 and 5 il 
cate that Snyder’s procedure yields results in fairly good agreement with t 
frequency data obtained from the observed streamflow record for the Rapp 
hannock and the synthesized floods for Goose Creek and Little Falls Brane 
However, the fact that in all three comparisons Snyder’s procedure yields 
highest 50- and 100-year floods, whereas his values for other recurrence 
intervals may be lower or higher, suggests a possible bias in his method. 

Incidentally, the writer was unable to determine how the author obtainec 
the Gumbel values he listed for Goose Creek and Little Falls Branch. The 
values listed above in Table 5 can be expected to be significantly higher thi 
the author’s for recurrence intervals of 10 years and under because they h 
been adjusted for a partial duration series. However, the differences are 
such that they cannot be explained by that alone. In both cases, the author’ 
computations appear to be at fault. In the case of Goose Creek, two annual 
peak discharges of 45,000 cfs and one of 32,800 cfs during the 26-year per: 
of record indicate that the author’s values are much too low. In the case o 
Little Falls Branch, the writer used an additional year of record but that 
value was the minimum annual peak discharge of record. Consequently, 
there should be no appreciable difference between the frequency values fre 
the 12-year record and those from the 13-year. However, the values pres 
ed by the author in his Table 2B are again appreciably lower than the writ 
(Table 5). Possible explanations are that the author obtained his Gumbel ¢ 
by fitting a line visually instead of by computation and/or that he did not et 
sider published and unpublished revisions of published annual peak discha1 

Perhaps the weakest link in the development of Snyder’s procedure is tl 
rainfall-runoff relation, which is based on monthly precipitation and runof 
data. He assumes that such a relation tends to represent storm situations 
advances questionable reasons, Nevertheless, he subjectively drew his cu 
to yield “slightly higher” runoff than would be indicated by a mean curve. 
Why? How much higher? One is inclined to conclude that the curve had t 
be drawn to yield some pre-determined results. 

The use of a single curve to represent even average storm rainfall-7al 
situations is not sound. For example, his relation (Fig. 3, App. 1) indicate 
about 48 per cent runoff for 5 inches of rainfall. Presumably, it makes no 
difference whether the 5-inch rainfall falls in one hour or 24 hours. This 
of course, contrary to fact and therefore unacceptable. A duration, or tim 
of concentration, parameter must be included in the relation. 

In order to determine how well Snyder’s relation fitted individual basin 
relations based on monthly precipitation-runoff data were derived for five 
the basins listed by him. Fig. A shows the relations for the smallest and’ 
largest natural basins, namely, Little Falls Branch and Rappahannock Riv 
with Snyder’s relation superimposed. In both cases, as in the other cases 
shown, a seasonal effect is noted. In the case of Little Falls Branch, Snyc 
relation appears to yield considerably higher than a mean curve for all th 
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MONTHLY PRECIPITATION ( In.) 
Figure A. MONTHLY RAINFALL-RUNOFF RELATIONS 


. would indicate. Fig. A suggests that Snyder’s mean relation does not 
serly represent mean monthly precipitation-runoff relations, let alone 

m rainfall-runoff conditions, for some individual basins in the area of 
rest. In fact, the indications are that there is a very small chance that 
relation would approximate average storm-runoff conditions on any basin. 
nother weak step in the development of Snyder’s procedure is the point- 
1 relation. He assumed a random distribution of storm centers over a 
ular basin and obtained two sets of reduction factors, one for summer 
ms and one for winter storms. The former was assumed to apply to 
-entration times of 4 hours and less andthe latter to 6 hours and more. 

e frequency is involved, it would have been much better to have used the 
t-area relation presented in U. S. Weather Bureau Technical Paper 29, 3 
hich the author refers. That relation is based on a large amount of ob- 
ed data of comparable frequencies from unusually dense precipitation 
orks and is therefore most applicable to Snyder’s problem. 

sS implied above, considerable subjectivity may be involved in the selection 
€ proper point rainfall-frequency data for the computation of flood 
uencies for a particular basin. Subjectivity may also be involved in the 
ction of a representative value of the friction factor, n, and in estimating 
ercentage of area that is impervious or drained by sewers and the 
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percentage of natural drainage channels eliminated. A difference of .01 in 
value of n may sometimes lead to about a 10-per cent difference in the con 
puted flood-frequency values. Errors of equal magnitude may result from 
incorrect evaluations of other factors used in determining Te. 

In summary, the chief criticisms of Snyder’s procedure apply to its de- 
velopment and not to the procedure itself. Some subjectivity was apparent 
involved in the development. This applies especially to the derivation of th 
rainfall-runoff relation. Furthermore, that relation and the point-area re- 
lation are oversimplified. Application of the procedure would also require 
some subjectivity in the selection of rainfall-frequency data, friction facto: 
etc., to be used. In view of these shortcomings, the writer cannot agree wi 
the author that his procedure would yield flood-frequency data of greater r 
liability than could be obtained by synthesizing flood records. (1) However, 
Snyder’s approach appears to be basically sound, and can be modified to yi 
results of much greater reliability than could be expected from it in its 
present form. 
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| DIVINING RODS VERSUS HYDROLOGIC DATA AND RESEARCH®@ 


| Discussion by Ven Te Chow 


VEN TE CHOW,! A. M. ASCE.—The author has presented an extremely 
mulating paper on the history, philosophy, and the general situation of the 
rologic data collection and treatment at the present time. This presen- 
on is timely and valuable because hydrology as either a science or a tech- 
ogy has come of age. Owing to a general negligence of the problem, the 

d for an adequate basic-data program in hydrology should be greatly 
ghasized and must be made known to all people concerned. 

The use of diving rods in the discussion of the problem in question is 
resting and appropriate. In fact, the philosophy of divining rods is not 
jue in the development of hydrologic problems, while its essence is also 

2 in many other scientific problems. The acceptance of witching philosophy 
the majority of the society constitutes indeed a stumbling-block to the 
gress of technology and to the development of any scientific program. 
According to a recent survey made in the study at Harvard University, a 
at number of approximately 25,000 diviners are known in existence, 

rently plying their trade in the U. S. A.2 The conclusion of this study indi- 
2s that “waterwitching will persist in this country as long as individuals 

st choose well sites under conditions where they have few external guides 
0 which of several possible sites will minimize their chances of paying 
vily for the water they need.” A modern basic-data program for collect- 
and interpreting hydrologic information should supply sufficient “external 
les” to those who are in need and therefore will end any attempt of a 

stical decision-maker. 

[The gap between the data collection and the data processing and interpreta- 
is apparent from the author’s statement: “About 20,000 pages of water 

| are published annually. Only a few hundred pages are devoted to inter- 
tive and research results.” It is economically unjustifiable that much ef- 
and hence money have been spent on data collection while relatively little 
rpretive results are produced. Although many agencies are engaging in 
collection and publication of hydrologic data, there is seemly lacking ina 
serted effort or unified attempt to make the information readily available 
researcher in the field. Many times the researcher has to spend much 
rgy in collecting the information from the information collected by differ- 
agencies. Therefore, one thing that a basic-data program should contain 
unified hydrologic-data agency, if proved to be feasible, should do is to 


Proc. Paper 1809, October, 1958, by W. B. Langbein. 

>rof. of Hydr. Eng., Univ. of Mlinois, Urbana, Il. 

tay Hyman and Evon Z. Vogt, 1958. Some Facts and Theories on Water- 
Vitching in the United States: GeoTimes II (9), March 1958, pp. 6-7 and 
9. F 
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classify and collect the information from the individual collections mad 
all public and private agencies. Such classification and collection can be 
easily recorded on IBM cards and an index to the information should be is 
from time to time. The U. S. Geological Survey is currently publishing i 
es of surface-water records in the United States. Such indexes should be 
panded and extended to cover all kinds of hydrologic data by a sort of cent 
ized agency. 
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WAVE FORCES ON SUBMERGED STRUCTURES2 


; Discussion by Turgut Sarpkaya 


CURGUT SARPKAYA, | A. M. ASCE.—The timeliness of this paper lies in 
fact that while it provides design data much needed for the application of 
lytical methods in the determination of forces on some types of submerged 
ictures, it shows the method of handling the elusive problem of determin- 
the drag in unsteady flow by means of an admittedly somewhat unrealistic 
quite simple analysis. 

[his writer’s comments and criticisms will be on testing facilities, test 
cedure, and on the basic concept of virtual mass and drag in unsteady 

la 


|. Under the heading “Test procedure”, the authors state that “For 

tain wave frequencies, standing transverse waves were developed which 
juced a variation in wave height from one side of the model to the other 
nuch as 40 per cent. For such conditions, a single measurement of wave 
tht at the center of the model provided a value which differed as much as 
per cent from the average effective wave height. This phenomenon was 
sed by the establishment of a resonant motion in the channel and tank, and 
nethod was found to eliminate it.” Basically the problem is that of the 
Sversal instability of oscillatory gravity waves in open channels. Suquet 
Wallet(1) and ites th sriscoolcnant reported similar instabilities. 

n(3) has demonstrated experimentally that when the frequency of the 
llations of the cylindric plunger exceeds a certain well defined critical 

e, transverse waves are also set up resulting in a three dimensional 
fitational oscillation. In addition to the resonant motion set up in the 

el basin as described by the authors, the presence of vertical walls of the 
nel and most important of all the particular wave generating mechanism 
elieved to be responsible for the transversal instability of waves. This 

ef is supported by the fact that a series of plates attached perpendicularly 
1e plunger at regular intervals have prevented or eliminated the trans- 

al standing waves. (1 

[the authors have employed such anti-standing wave plates on the plunger, 
> filters immediately at the downstream side of the plunger, and artificial 
2 absorbers at the end of the test channel, they might have prevented the 
lation of standing transverse waves. 

. This writer believes that the authors are well justified in using Airy’s 
tion for oscillatory gravity waves, rather than those of Stokes or Struik, 
he reason that the experimental results reported(4,5) are in better 


Sd 

roc. Paper 1833, November, 1958, by Ernest F. Brater, John S. McNown, 
nd Leslie D. Stair. 
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agreement with Airy’s equation as far as the velocity of waves are con- 
cerned. However, it should be noted that the wave profile is described bett 
by Stokes’ or Struik’s equations. The deviation of a sinusoidal curve from 
the actual profile might slightly effect the AH value. 

3. Computation of oscillatory wave forces or the sum of the drag and th 
inertial force on the submerged structures under consideration is based up 
the assumption that the flow pattern is not disturbed by the presence of the 
model. For the flat plates the values of Cg and Cm were selected to give # 
best correspondence with the measured values of phase angles and the mag 
tudes of the maximum forces. For barge models the values of phase angle 
and the magnitudes of the maximum forces. For barge models the values ¢ 
Cm were computed using the measured values of inertial forces. Hence, fc 
the barges there is'no way of knowing how approximate the values of Cy a1 
and how much justification there is for the simplified analysis. One can on 
say that Cy, is assumed to take into account whatever discrepancies may 
arise between a more realistic analysis and a simplified one. For rectang 
lar barge models the authors found a Cy value which varied from 1.31 to 
1.76. At first these values seem to compare favorably with the theoretical 
value of C,;, obtained from a two-dimensional flow analysis ignoring the ef 
fects of the time dependent wake surrounding the body. If, however, one pi 
pares a plot of the rate of change of Cy) with respect to the ratio (a/b), it i 
easy to observe that for the barge model tested the value of C,, does not v: 
more than a few per cent if only the width of the model is increased by 509% 
i.e., from 10 in., to 15 inches. Hence, a wake that would form only behind 
plate would not have any appreciable effect on the value of Cm. Since, how 
ever, the sharp edges lead to separation for very small relative motions, 1 
sulting vortices grow into a wake and extend the disturbed flow region also 
laterally as well as transversely. Therefore, in comparing the experimen’ 
Cy values with the ones obtained from a two dimensional flow analysis, on 
should take into consideration the equivalent potential flow model which is 
supposed to be made of a surface enclosing both the body and its time de- 
pendent wake.. 

In order to determine the effect of the three dimensional character of t] 
models on the value of C,,, plexiglass models (1/2 geometrical ratio) of th 
right rectangular prism and the slotted barges shown in the authors’ Fig. + 
were constructed in the hydrodynamics laboratory of the Engineering Me- 
chanics department of the University of Nebraska, and the C values (added 
mass/displaced mass) were determined by means of vibratory motion. Fe 
full rectangular barge the value was of 0.35, and for the slotted barge 0.65 
Corresponding C,, values are 1.35, and 1.65, respectively. In the tests at 
Nebraska, the velocities varied from 0.2 to 0.6 in/sec., the frequencies fr 
20 to 60 cycles per second, and the accelerations from 40 to 400 in/sec.se 
The above values of Cy, will be used later in computing the reduction in he 
zontal forces on slotted barges. 

In what follows, this writer would like to review the recent ideas conce 
ing virtual mass and then point out the parameters that cause marked char 
in both Cq and Cm. Although, considerable progress has been made on the 
experimental techniques and although the added mass coefficient has been 
termined for a wide range of shapes of bodies, no new physical meaning hé 
been attached to it, and the basic definition of virtual mass has remained ; 
“The quotient of the force required to produce the accelerations throughou 
the fluid divided by the acceleration of the body.” The usual derivation of 


DISCUSSION 69 


al mass is from kinetic energy which suffers from the inherent weakness 
it assumes a constant velocity for the confined solid body, whereas all 
yard manifestations of mass should be associated with acceleration. 
refore, the velocity must be allowed to vary and the fact should be recog- 
d that the added mass can yield momentum as well as energy. 
urprisingly enough, as early as 1888, Riecke(6) recognized the essential 
jthat for a sphere moving with constant velocity through an infinite 
scid fluid, the individual fluid particles which are pushed aside by the 
re in its forward motion do not return to their former positions. The 
is of the individual particles are not closed curves but of the shape simi- 
ro that shown in Fig. 1, of this paper. Hence, besides pushing the parti- 
aside temporarily in passing, the sphere also displaces the fluid parti- 
permanently in the direction of its motion. However, the importance of 
displacement and its relation to added mass was not recognized by 
‘ke. Recently, Darwin(7) has shown that this permanently displaced mass 
1e fluid enclosed between the initial and final positions of fluid particles, 


FIG.I 


‘R is cylinder's radius. Trajectories are shown 
‘for three particles. Broken lines show the in- 
‘itial and final positions of particles before 
and after the passage of the cylinder 
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is in fact the added mass itself. It can be shown, in general, that the motic 
of a body through an inviscid fluid media is always accompanied by a fluid- 
mass transport and that this mass is the added mass which unveils itself 
only if the body is accelerated. 

The above explanation, although adding considerably to our understandir 
of the added mass, does not take into account the effect of the wake- or the 
cavity-induced mass. Furthermore, for an accelerating or decelerating be 
the induced mass varies with the instantaneous shape and volume of the wa 
or cavity as well as with their rates of change, i.e., the instantaneous magi 
tude of added mass values in transient conditions depends upon the time 
history of motion. Suggestions have been made to approximate transient 
motions by parts of a sine curve and thus to determine a frequency paramée 
ter. This procedure is obviously questionable. It is clear from the above 
discussion that the variation of the wake geometry with time causes marke 
changes in both Cg and Cy). Because both coefficients depend upon the sta 
of development of the wake, recently McNown, 8) McNown and Keulegan 9 
lated them to each other quite successfully. 

It should be noted, furthermore, that in an oscillatory gravity wave the 
elliptical orbits are open and there is a continuous mass transport in the ¢ 
rection of wave propagation near the surface, and in the opposite direction 
near the bottom. Therefore, the value of C,, computed from experimental 
data also involves the effect of the inertia of wave mass transport. 

4. In order to obtain the pressure difference on the two sides of the 
volume of fluid which replaces the submerged structure, the authors could 
have used their Eq. (11), together with Eq. (4). If this is done for the case 
a rectangular barge centered at the quarter point of a wave, one obtains 


Li 
4 2, 
y a7 2 Ee OE ppd Be 27H cosh 27(4 + 2)/L 54, 22% gy 
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or performing the simple integration, 
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in which a, and ay represent respectively the lengths of horizontal axes of 
water particle orbits at a point z, and at a point on the still water level. _ 

If Eq. (2) is compared with the authors’ Eq. (16), it will be noticed thai 
second term in the parenthesis of the latter equation is missing in the for 
equation. The difference is basically due to the linearized averaging of th 
horizontal accelerations along the submerged structure. An exact analysi 
made by this writer of the difference in the pressures occurring under a ‘ 
crest and a trough using Struik’s equations(10) has shown that the right sit 
of the authors’ Eq. (18), contains more negative terms in the form of a se) 
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apidly decreasing absolute magnitudes. This analysis is not given here 
because it is irrelevant to this discussion and partly because it is 
> lengthy. 
ee authors have found that the reduction in displaced mass of the 
by introducing a slot is partially offset by, “the additional disturbance 
utable to the odd configuration of the slotted barge . . .”, and the values 
_ are greater than the ones for a full rectangular barge by 10 to 24 per 
i: values of H/L from 0.02 to 0.06. 
study, made by this writer, to determine the virtual mass of multiple 
es, lens shaped bodies, and of circular, square, and rectangular plates 
rious sizes amply supports only the results of the above observation, 
pt the explanation given for its cause. Since the results of this study 
ie reported later, only part of the results obtained for parallel square 
sand for sufficiently long (large ratio of length to width) parallel rec- 
lar plates will be given here. 
Figs. 2 and 3, curves with dotted lines give the ratio (added mass/ 
iced mass) for full right square prism and for unslotted right rectangu- 
‘ism, as a function of (s/a) for constant values of (t/a). Curves with full 
give the aforementioned ratio for the separated square and rectangular 
3 as a function of the same dimensionless ratios. Here, only the repre- 
tive curves, rather than the separate experimental points, are shown 
e sake of simplicity. The direction of vibratory motion is shown on the 
Ss. 
examination of the figures reveals that as the plates are separated 
each other but a very little distance, the added mass coefficient de- 
es slightly, i.e., the fluid between the plates behaves as if it were part 
Solid full parallelepiped. As the plates are set farther apart, the add- 
iss coefficient increases, but always remains smaller than the added 
coefficient of a single plate. Hence, when the distance between the 
is not too large in comparison with the width of the plate, there is a 
1 interaction of the flow patterns around the two plates. As the plates 
1e farther apart, this interaction becomes negligible and each plate be- 
independently. 
r the slotted barge model used by the authors t/a = 1.265, and s/a = 4. 
point corresponding to these parameters is located on Fig. 3, it be- 
3 immediately apparent that each leg of the slot behaves practically as 
ependent parallelepiped. Therefore, the reason for the partial offset of 
duction in displaced mass and the resulting increase in Cy, is not be- 
of the additional disturbance attributable to the odd configuration of the 
i barge, but on the contrary is because of the lack of interference of the 
des of the slot. For the forces in the horizontal direction, a weighted 
ye for Cy, can be determined easily. For the unslotted portion 
1.33, and for each leg Cm = 1.95, and the weighted average, although 
<imate for this shape of body, is 1.55. This value is in the neighborhood 
t has been found by the authors for various values of T, z/d, and H/L. 
ct that the length of the legs is not infinite and that the slotted barge 
is actually more complicated is, of course, not overlooked. But, this 
Mt alter the above reasoning. Using the Cy, values determined experi- 
lly by this writer for the two types of barges (Cy = 1.65 for the slotted 
‘and Cy = 1.35 for the rectangular barge) a reduction of only 9.5 per 
| the horizontal force is computed. 


March, 1959 


72 


ie ee aE wer Ce Ur 


WG “SiG wpe 


a 


4 


ssDW poaovidsip 
—s50u peppo 


DISCUSSION 


HY 3 


¢°9l4 


Ov gf OF ve Ze O'€ BZ O2 weZ SenOS Bs Siler eects Oe ee 20.0 


D ie) 
“| | | a scpees eel yelled htop} — 
Ce ee Se 
Pe mee =a : 
zi ee > Raat ee = 
— ios me 


714 


March, 1959 HY 3 


velocity (Eqs. (4) and (5)), and the equations for horizontal and vertical ac- 
celerations (Eqs. (6) and (7)), and hence, the tabulated values in Fig. 5, are 
not mathematically consistent as far as the rules of differentiation are con- 
cerned and as long as the meaning attached to theta is not changed in the 
meantime. 
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#98(PO2), 1599(PO2), 1600(PO2), 1601(P02), 1602(P02), 1603(HY2), 1604(EM2), 1605(SU1)¢,1606(SA2), 
HOT(SA2), 1608(SA2), 1609(SA2), 1610(SA2), 1611 (SA2), 1612(SA2), 1613(SA2), 1614(SA2)°, 1615(1R2)°, 1616 
T¥2)¢, 1617(SU1), 1618(PO2)°, 1619(EM2)°, 1620(CP1). 


7 1621(HW2), 1622(HW2), 1623(HW2), 1624(HW2), 1625(HW2), 1626(HW2), 1627(HW2), 1628(HW2), 1629 
'T3), 1630(ST3), 1631(ST3), 1632(ST3), 1633(ST3), 1634(ST3), 1635(ST3), 1636(ST3), 1637(ST3), 1638(ST3), 
38(WW3), 1640(WW3), 1641(WW3), 1642(WW3), 1643(WWS3), 1644(WW3), 1645(SM2), 1646(SM2), 1647 
M2), 1648(SM2), 1649(SM2), 1650(SM2), 1651(HW2), 1652(HW2)¢, 1653(WW3)°, 1654(SM2), 1655(SM2), 
5556 (ST3)°, 1657(SM2)°. 
fb: 1658(AT1), 1659(AT1), 1660(HY3), 1661(HY3), 1662(HY3), 1663(HY3), 1664(HY3), 1665(SA3), 1666 
»L2), 1667(PL2), 1668(PL2), 1669(AT1), 1670(P03), 1671(PO3), 1672(PO3), 1673(PL2), 1674(PL2), 1675 
703), 1676(P03), 1677(SA3), 1678(SA3), 1679(SA3), 1680(SA3), 1681(SA3), 1682(SA3), 1683(PO3), 1684 
MTY3), 1685(SA3), 1686(SA3), 1687(P.03), 1688(SA3)°, 1689(P.03)°, 1690(HY3)°, 1691(PL2)°. 


2) 
Fr: 1692(EM3), 1693(EM3), 1694(ST4), 1695(ST4), 1696(ST4), 1697(SU2), 1698(SU2), 1699(SU2), 1700(sU2), 
 01(SA4), 1702(SA4), 1703(SA4), 1704(SA4), 1705(SA4), 1706(EM3), 1707(ST4), 1708(ST4), 1709(ST4), 1710 
4), 1711(ST4), 1712(ST4), 1713(SU2), 1714(SA4), 1715(SA4), 1716(SU2), 1717(SA4), 1718(EM3), 1719 
\MB}, 1720(SU2), 1721(ST4)°, 1722(ST4), 1723(ST4), 1724(EM3)°. 
us 3), 1732(SM3), 173 
. 1725(H¥4), 1726(HY4), 1727(SM3), 1728(SM3), 1729(SM3), 1730(SM3), 1731(SM3), 
ta 1734(PO4), 1735(PO4), 1736(PO4), 1737(PO4), 1738(PO4), 1739(PO4), 1740(PO4), 1741(PO4), 1742 
QA), 1743(PO4), 1744(P04), 1745(PO4), 1746(PO4), 1747(PO4), 1748(PO4), 1749(PO4). 


: 1750(0R3), 1751(0R3), 1752(1R3), 1753(0R3), 1754(1R3), 1755(ST5), 1756(ST5), 1757(ST5), 
: fare) See em 3 1761(STS), 1762(STS), 1763(ST5), 1764(STS), 1765(Ww4), 1766(WW4), 
-en(ww4), 1768(WwWw4), 1769(WW4), 1770(WW4), 1771(WW4), 1772(WW4), 1773(WWw4), 1774(IR3), 1775 
1776(SA5), 1777(A5), 1778(SA5), 1779(SA5), 1780(SA5), 1781(WW4), 1782(SA5), 1783(SA5), 1784 
fe, 1785(wwA)®, 1786(SA5)°, 1787(STS)°, 1788(1R3), 1789(WW4). aie 
BER: 1791(EM4), 1792(6M4), 1793(EM4), 1794(2M4), 1795(HW3), 1796(HW3), 1797(HW3), 
ecw), ree ence oe 1BOREW 3) ecadawa), repa(nw), 1B0503), 1806 
r5), 180" 1808(HY5), 1809(HY5), 1810(HY5), 1 ‘ ; : : 
aan fseen. ee 1819(ST6), 1820(ST6), 1621 (S78), 1822(EMA), 1823(PO5), 1824 
4), 1825(SM4), 1826(SM4), 1827(STS)°, 1828(SM4)°, 1829(HW3)°, 1830(PO5)°, 1831(EM4)°, 1832(HY5) 


=MBE 2 : 1834(HY6), 1835(SA6), 1836(ST7), 1837(ST7), 1838(ST7), 1839(ST7), 1840(ST7), 
343 rn, “1e4a@US), zeAs GUS, 1844(SU3), 1845(SU3), 1846(SU3), 1847(SA6), 1848(6A8), 149(6A6), 
350(SA6), 1851(SA6), 1852(SA6), 1853(SA6), 1854(ST7), 1855(SA6)°, 1856(HY6)°, (sT7)°, 


MBER: 1859(HY7), 1860(1R4), 1861(IR4), 1862(1R4), 1863(SM5), 1864(SMS), 1865(ST8), 1866(STS), 1867 

eeierts pie. hen 1871(PP1), 1872(PP1), 1873(WW5), 1874(WW5), 1875(WWS), 1876 
NWS), 1877(CP2), 1878(ST8), 1879(ST8), 1880(HY7)°, 1881(SM5)°, 1882(ST8)°, 1883(PP1)°, 1884(Wws)°, 
#885(CP2)¢, 1886(PO6), 1887(PO6), 1888(PO6), 1889(PO6), 1890(H¥7), 1891(PP1). 


VOLUME 85 (1959) 


JARY: , 1894(EM1), 1895(BM1), 1896(EM1), 1897(EM1), 1898(EM1), 1899(HW1), 
30 LN Sater te Ubeastaldeg eae 1904(HY1), 1905(PL1), 1906(PL1), 1907(PL1), 1908(PL1), 
09(ST1), 1910(ST1), 1911(ST1), 1912(ST1), 1913(ST1), 1914(6T1), 1915(8T1), 1916(AT1)°, aM 
(HW1)°, 1919(HY1)°, 1920(PL1)°, 1921(SA1)°, 1922(ST1), 1923(EM1), 1924(HW1), 1925(HW1), 
27(HW1), 1928(HW1), 1929(SA1), 1930(SA1), 1931(SA1), 1932(SA1). 
RY: 1933(HY2), 1934(HY2), 1935(HY2), 1936(SM1), 1937(SM1), 1938(ST2), 1939(ST2), 1940(ST2), 
2), 1942(ST2), 1943(ST2), 1944(ST2), 1945(HY2), 1946(PO1), 1947(PO1), 1948(PO1), faAsPou), 
TY2)¢, 1951(SM1)¢, 1952(ST2)°, 1953(PO1)S, 1954(COL), 1955(COL), 1956(CO1), 1957(CO1), 19 
959(COl). 
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